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The string-net approach by Levin and Wen, and the local unitary transformation approach by 
Chen, Gu, and Wen, provide ways to systematically describe topological orders in 2D boson sys- 
tems. The two approaches reveal that the mathematical framework for topological order is closely 
related to (projective) tensor category theory. In this paper, we generalize those systematic de- 
scriptions of topological orders to 2D fermion systems. We find that a large class of topologi- 
cal orders in 2D fermion systems can be systematically characterized by the following set of data 
(Nijk, Nfj k , Fjj^'^sC i Ofa^ 13 ) that satisfy a set of non-linear algebraic equations. The exactly solu- 
ble Hamiltonians can be constructed from the above data on any lattices to realize the corresponding 
topological orders. Since topological orders in 2D boson systems described by the string-net ap- 
proach have a deep relationship with the Turaev-Viro states summation invariance of 3-manifolds, 
it is natural to expect our classification of topological orders in 2D fermionic systems will generalize 
the famous Turaev-Viro invariance to its fermionic version. When N( jk = 0,' our results provide a 
classification of bosonic topological orders that generalizes the string-net classification. 



3 



I 

-a 

G 

O 
O 



> 

in 



o 
o 



X 



CONTENTS 



I. Introduction 



II. Local fermion systems 3 

A. Local bosonic operators and bosonic states in a 
local boson model 3 

B. Local fermionic operators and fermionic states in 

a local fermion model 3 

III. Fermionic local unitary transformation and 
topological phases of fermion systems 3 

IV. Fermionic local unitary transformation and wave 
function renormalization 4 



V. Wave function renormalization and a classification 
of fermionic topological orders 

A. Quantum state on a graph 

B. The structure of a fixed-point wave function 

C. The first type of wave function renormalization 

D. The second type of wave function 
renormalization 

E. The fixed-point wave functions from the 
fixed-point gfLU transformations 

F. Summary of fixed-point gfLU transformations 

VI. Generalized Turaev-Viro invariance and its graphic 
representation 

A. The amplitude of tetrahedra graph and 
fermionic symmetric Q symbol 

B. Graphic representation for Q symbol 

1. Graphic representation in the bosonic case 

2. Graphic representation for the full Q symbol 

C. The 2^3/3^2 move for the Q symbol 

1. The orthogonality condition 

2. Graphic representation for the symbol 

3. The 2 -> 3/3 -> 2 move 

D. The 2 -> 0/0 -> 2 move for the Q symbol 
1. The associativity condition of A 1 



5 
5 
5 
6 

11 

12 
13 



14 

14 
15 
15 
15 
16 
16 
16 
16 
17 
17 



2. The 2^0/0^2 move 18 

E. Topologically invariant partition function 19 

F. Simplified fixed-point conditions 19 

VII. Categorical framework 20 

A. Projective super tensor category 20 

B. Super tensor category from super quantum 
groups 20 

VIII. Simple solutions of the fixed-point conditions 20 

A. Unimportant phase factors in the solutions 21 

B. A Z3 bosonic fixed-point 21 

C. A Z3 fermionic fixed-point 22 

D. A class of "trivial" fermionic fixed points 25 

IX. Summary 25 

A. The fermionic structure of support space 25 

B. Ideal Hamiltonian for the fixed-point states 25 

C. Putting fermions on both vertices and edges 30 

1. Quantum state on a graph 30 

2. The first type of wave function renormalization 30 

3. The second type of wave function 
renormalization 32 

4. Fixed-point conditions on the F-tensor and the 
O-tesnor 32 



References 



33 



I. INTRODUCTION 

Understanding phases of matter is one of the central 
problems in condensed matter physics. Landau symme- 
try breaking theory, 1 ' 2 as a systematic theory of phases 
and phase transitions, becomes a corner stone of con- 
densed matter theory. However, at zero temperature, 
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the symmetry breaking states described by local order 
parameters are basically direct product states. It is hard 
to believe that various direct product states can describe 
all possible quantum phases of matter. 

Based on adiabatic evolution, one can show that 
gapped quantum phases at zero temperature correspond 
to the equivalence classes of local unitary (LU) trans- 
formations generated by finite-time evolutions of local 
hermitian operators H(t): 3 ~ 6 

|$) ~ |$'> iff |*') =Te i / dT ^ T )|$). (1) 

It turns out that there are many gapped quantum states 
that cannot be transformed into direct product states 
through LU transformations. Those states are said to 
have a long range entanglement. Thus, the equivalence 
classes of LU transformations, and hence the quantum 
phases of matter, are much richer than direct prod- 
uct states and much richer than what the symmetry 
breaking theory can describe. Different patterns of long 
range entanglement correspond different quantum phases 
that are beyond the symmetry-breaking/order-parameter 
description 7 and direct-product-state description. The 
patterns of long range entanglement really correspond 
to the topological orders 8 ' 9 that describe the new kind 
of orders in quantum spin liquids and quantum Hall 
states. 10 " 18 

In absence of translation symmetry, the above LU 
transformation can be expressed as a quantum circuit, 
which corresponds to a discretized LU transformation. 
The discretized LU transformation is more convenient 
to use. The gapped quantum phases can be more effec- 
tively studied and even classified through the discretized 
LU transformations. 6 ' 19 ~ 21 

After discovering more and more kinds of topologi- 
cal orders, it becomes important to gain a deeper un- 
derstanding of topological order under a certain math- 
ematical framework. We know that symmetry break- 
ing orders can be understood systematically under the 
mathematical framework of group theory. Can topologi- 
cal orders be also understood under some mathematical 
framework? From the systematic construction of topo- 
logically ordered states based on string-nets 19 and the 
systematic description of non-Abelian satistics 22 , it ap- 
pears that tensor category theory may provide the under- 
lying mathematical framework for topological orders. 23 

However, the string-net and the LU transformation 
approaches 6 ' 19 only provide a systematic understanding 
for topological orders in qubit systems (ie quantum spin 
systems or local boson systems). Fermion systems can 
also have non trivial topological orders. In this paper, 
we will introduce a systematic theory for topological or- 
ders in interacting fermion systems (with interacting bo- 
son systems as special cases). Our approach is based 
on the LU transformations generated by local hermitian 
operators that contain even number of fermion opera- 
tors. It allows us to classify and construct a large class 
of topological orders in fermion systems. The mathe- 
matical framework developed here may be related to the 



theory of enriched categories, 2 which can be viewed as a 
generalization of the standard tensor category theory 25 ' 26 
To gain a systematic understanding of topological or- 
der in fermion systems, we first need a way to la- 
bel those fermionic topological orders. In this pa- 
per, we show that a large class of fermionic topolog- 
ical orders (which include bosonic topological orders 
as special cases) can be labeled by a set of tensors: 

(N ijk ,N{ jk ,F;i^C,Oi k a^) : Certainly, not every set 
of tensors corresponds to a valid fermionic topological or- 
der. We show that only the tensors that satisfy a set of 
non-linear equations correspond to valid fermionic topo- 
logical orders. The set of non-linear equations obtained 
here is a generalization of the non-linear equations (such 
as the pentagon identity) in a tensor category theory. So 
our approach is a generalization of tensor category the- 
ory and the string-net approach for bosonic topological 
orders. We like to point out that the framework devel- 
oped here not only lead to a classification of fermionic 
topological orders, it also lead to a more general classi- 
fication of bosonic topological orders than the string-net 
and the related approaches. 6 ' 19 

From a set of the data (N ijk ,N( jk ,Fj^ a , Oft""), 
we can obtain an Hamiltonian as a sum of projectors. We 
believe that the Hamiltonian is unfrustrated. Its zero- 
energy ground state realizes the fermionic/bosonic topo- 
logical order described by the data. For a subclass of the 
topological order, the data satisfies a set of non-linear 
equation which are not projective. The ideal Hamiltoni- 
ans for such a subclass of topological orders are sums of 
commuting projectors, and are thus exactly soluble. 

In this paper, we will concentrate on bosonic and 
fermionic topological orders in 2D interacting systems. 
We like to mention that recently some significant pro- 
gresses have been made for bosonic and fermionic (sym- 
metry protected) topological orders in ID interacting sys- 
tems: a complete classification is obtained, 27-29 which 
generalizes some partial results obtained earlier. 20 ' 30 ~ 34 

In section II, we give a careful discussion on what 
is a local fermion system. In section III, we introduce 
fermionic local unitary transformations. We then use 
fermionic local unitary transformations to define quan- 
tum phases for local fermion systems and fermionic topo- 
logical orders. In section IV, we use fermionic local uni- 
tary transformations to define a wave function renormal- 
ization flow for local fermion systems. In section V, we 
discuss the fixed points of the wave function renormal- 
ization flow, and use those fixed points to classify a large 
class of fermionic (and bosonic) topological orders. In 
section VI, we present a graphic representation of some 
of our results. In section VII, we comment on the rela- 
tion to categorical framework. In section VIII, we give 
a few simple examples. In appendix A, we discuss the 
fermionic structure of the support space. In appendix 
B, we derive the ideal Hamiltonian from the data that 
characterized the fermionic/bosonic topological orders. 
In appendix C, we describe an even more general classi- 
fication of fermionic topological orders. 



3 



II. LOCAL FERMION SYSTEMS 

Local boson systems (ie local qubit systems) and local 
fermion systems have some fundamental differences. To 
reveal those differences, in this section, we are going to 
define local fermion systems carefully. To contrast local 
fermion systems with local boson systems, let us first 
review the definition of local boson systems. 



A. Local bosonic operators and bosonic states in a 
local boson model 

A local boson quantum model is defined through its 
Hilbert space V and its local boson Hamiltonian 77. The 
Hilbert space V of a local boson quantum model has a 
structure V — <g)Vi where Vi is the local Hilbert space 
on the site i. A local bosonic operator is defined as an 
operator that act within the local Hilbert space Vi, or 
as a finite product of local bosonic operators acting on 
nearby sites. A local boson Hamiltonian H is a sum 
of local bosonic operators. The ground state of a local 
boson Hamiltonian H is called a bosonic state. 



B. Local fermionic operators and fermionic states 
in a local fermion model 

Now, let us try to define local fermion systems. A 
local fermion quantum model is also defined through its 
Hilbert space V and its local fermion Hamiltonian Hf. 
First let us introduce fermion operator cf at site i as 
operators that satisfy the anticommutation relation 

c ? c i = - c i c ?> c i ( c f) ! = -( c i) tc ?> ( 2 ) 

for all i ^ j and all values of the a, (3 indices. We also 
say c" acts on the site i. The fermion Hilbert space V is 
the space generated by the fermion operators and their 
hermitian conjugate: 

V = {[<£(c?)t...]|0>}. (3) 

Due to the anticommutation relation (2), V has a form 
V — <E)Vi where Vi is the local Hilbert space on the site i. 
We see that the total Hilbert space of a fermion system 
has the same structure as a local boson system. 

Using the Hilbert space V = ®V£, an explicit represen- 
tation of the fermion operator cf can be obtained. First, 
each local Hilbert space can be splitted as Vi = V? © V^ . 
We also choose an ordering of the site label i. Then c" 
has the following matrix representation 

c? = cy JjE|, 

C ? = \Bf A Q ) ' E * = ( -I}) ' (4) 



where 7? is the identity matrix acting in the space V® 
and l\ is the identity matrix acting in the space ■ 
The matrix Cf maps a state in VP to a state in U/ , and 
vice versa. We note that 

C?£? = -EfC? . (5) 

We see that a fermion operator is not a local bosonic 
operator. The product of an odd number of fermion op- 
erators and any number of local bosonic operators on 
nearby sites is called a local fermionic operator. 

Let us write the eigenvalue of S? as (— ) Si . The states 
in VP have Si = and are called bosonic states. The 
states in V^ have Si = 1 and are called fermionic states. 
We can view Si as the fermion number on site i. 

A local fermion Hamiltonian Hf is a sum of terms: 
77 = ^ZpOp, where J2p sums over a set of regions. 
Each term Op is a product of an even number of lo- 
cal fermionic operators and any number of local bosonic 
operators on a finite region P. Such kind of terms is 
called pseudo-local bosonic operator acting on the region. 
In other words, a local fermion Hamiltonian is a sum of 
pseudo-local bosonic operators. The ground state of a 
local fermion Hamiltonian 77/ is called a fermionic state. 

Note that, beyond ID, a pseudo-local bosonic opera- 
tor is in general not a local bosonic operator. So a lo- 
cal fermion Hamiltonian 77/ (beyond ID) in general is 
not a local boson Hamiltonian defined in the last sub- 
section. In this sense, a local boson system and a local 
fermion system are fundamentally different despite they 
have the same Hilbert space. When viewed as a boson 
system, a local fermion Hamiltonian corresponds to a 
non-local boson Hamiltonian (beyond ID). Thus classi- 
fying the quantum phases of local fermion systems corre- 
sponds to classifying the quantum phases of a particular 
kind of non-local boson systems. 



III. FERMIONIC LOCAL UNITARY 
TRANSFORMATION AND TOPOLOGICAL 
PHASES OF FERMION SYSTEMS 

Similar to the local boson systems, the finite-time evo- 
lution generated by a local fermion Hamiltonian defines 
an equivalence relation between gapped fermionic states: 

mi)) ~ \m) iff \m) = rw ^^)]|^(o)) (6) 

where T is the path-ordering operator and 77(g) = 
^2iOi(g) is a local fermion Hamiltonian (ie Oi(g) is a 
pseudo-local bosonic operator which is a product of even 
local fermionic operators). We will call T[e'^o d s H f(9)^ 
a fermion local unitary (fLU) evolution. We believe that 
the equivalence classes of such an equivalence relation is 
the universality classes of the gapped quantum phases of 
fermion systems. 

The finite-time fLU evolutions introduced here is 
closely related to fermion quantum circuits with finite 
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depth. To define fermion quantum circuits, let us in- 
troduce piece-wise fermion local unitary operators. A 
piece-wise fermion local unitary operator has a form 
Up W i = Yii e lH fW, where {Hf(i)} is a Hermitian op- 
erator which is a pseudo-local bosonic operator that acts 
on a region labeled by i. Note that regions labeled by 
different i's are not overlapping. Ui = e, lH fW is called a 
fermion unitary operator. The size of each region is less 
than some finite number I. The unitary operator U pw i 
defined in this way is called a fermion piece-wise local 
unitary operator with range I. A fermion quantum circuit 
with depth M is given by the product of M fermion piece- 
wise local unitary operators: U^ rc = U^U^ ■ ■ ■ . 
We believe that finite time fLU evolution can be simu- 
lated with a constant depth fermion quantum circuit and 
vice versa. Therefore, the equivalence relation eqn. (6) 
can be equivalently stated in terms of constant depth 
fermion quantum circuits: 



i^(i)) ~ \m) iff \m) = u™ c \m) 



(7) 



where M is a constant independent of system size. Be- 
cause of their equivalence, we will use the term "fermion 
Local Unitary Transformation" to refer to both fermion 
local unitary evolution and constant depth fermion quan- 
tum circuit in general. 

Just like boson systems, the equivalence classes of 
fermionic local unitary transformations correspond to 
the universality classes that define phases of matter. 
Since here we do not include any symmetry, the equiva- 
lence classes actually correspond to topologically ordered 
phases. Such topologically ordered phases will be called 
fermionic topologically ordered phases. 



IV. FERMIONIC LOCAL UNITARY 
TRANSFORMATION AND WAVE FUNCTION 
RENORMALIZATION 

After defining the fermionic topological order as the 
equivalence classes of many-body wave functions under 
fLU transformations, we like to use the fLU transfor- 
mations, or more precisely the generalized fermion local 
unitary (gfLU) transformation, to define a wave function 
renormalization procedure. The wave function renormal- 
ization can remove the non-universal short-range entan- 
glement and make generic complicated wave functions to 
flow to some very simple fixed-point wave functions. The 
simple fixed-point wave functions can help us to classify 
fermionic topological orders. 

Let us first define the gfLU transformation U g more 
carefully. Consider a state \ip). Let pa be the entangle- 
ment density matrix of \ip) in region A. pa may act in a 
subspace of the Hilbert space in region A. The subspace 
is called the support space Va of region A (see Fig. 1(a)). 
Let \(j>i) be a basis of this support space Va, and be a 
basis of the full Hilbert space Va of region A. The gfLU 
transformation U g is the projection from the full Hilbert 
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FIG. 1. (Color online) (a) A fermion local unitary (gfLU) 
transformation U g acts in region A of a fermionic state |i/>) 
are formed by bosonic and fermionic operators in the region 
A. U g always contain even numbers of fermionic operators, 
(b) U^Ug = P is a projector, whose action does not change 
the state lib). 



space Va to the support space Va- So up to some unitary 
transformations, U g is a hermitian projection operator: 



Ug = U 1 PgU 2 , 

UlUi = 1, 



P A — p 
u\u 2 = 1. 



pt _ p 



(8) 



The matrix elements of U g are given by (<j)a\4>i). We will 
call such a gfLU transformation a primitive gfLU trans- 
formation. A generic gfLU transformation is a product of 
several primitive gfLU transformations which may con- 
tain several hermitian projectors and unitary transfor- 
mations, for example, U g — UxPgUzP'Uz. 

To understand the fermionic structure of U g , we note 
that the support space Va has a structure Va = V% © V\ 
(see appendix A), where V^ has even numbers of fermions 
and V\ has odd numbers of fermions. This means that U g 
contains only even numbers of fermionic operators (ie U g 
is a pseudo-local bosonic operator). 

We also regard the inverse of U g , t/J, as a gfLU trans- 
formation. A fLU transformation is viewed as a special 
case of gfLU transformations where the degrees of free- 
dom are not changed. Clearly U^Ug = P and U g W g = P' 
are two hermitian projectors. The action of P does not 
change the state \4>) (see Fig. 1(b)). Thus despite the 
degrees of freedom can be reduced under the gfLU trans- 
formations, no quantum information of the state is 
lost under the gfLU transformations. 

We note that the gfLU transformations can map one 
wave function to another wave function with fewer de- 
grees of freedom. Thus it can be viewed as a wave func- 
tion renormalization group flow. If the wave function 
renormalization leads to fixed-point wave functions, then 
those fixed-point wave functions can be much simpler, 
which can provide an efficient or even one-to-one label- 
ing scheme of fermionic topological orders. 
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FIG. 2. (Color online) A quantum state on a vertex-labeled 
graph G. The four vertices are labeled by a, f3, \, and 5. 
For example, a = 1, ,3 = 2, \ = 3, 5 = 4, correspond to 
one vertex labeling scheme. States on each edge are labeled 
by I = 0, N . If the three edges of a vertex are in the states 
i, j, and k respectively, then the vertex has Nijk states. The 
states on the a-vertex are labeled by a — l,...,Nijk- Note 
the orientation of the edges are point towards to vertex. Also 
note that i — > j — s> k runs anti-clockwise. 



FIG. 3. The mapping i —, 
the orientation of the edge 



corresponds to the reverse of 



V. WAVE FUNCTION RENORMALIZATION 
AND A CLASSIFICATION OF FERMIONIC 
TOPOLOGICAL ORDERS 



As an application of the above fermionic wave func- 
tion renormalization, in this section, we will study the 
structure of fixed-point wave functions under the wave 
function renormalization. This will lead to a classifica- 
tion of fermionic topological orders. 



A. Quantum state on a graph 

Since the wave function renormalization may change 
the lattice structure, we will consider quantum state de- 
fined on a generic trivalent graph G: Each edge has N + 1 
states, labeled by i = 0, ...,7V (see Fig. 2). We assume 
that the index i on the edge admits a one-to-one mapping 
*: i —> i* that satisfies (i*)* = i. As a result, the edges of 
the graph are oriented. The mapping i — > i* corresponds 
to the reverse of the orientation of the edge (see Fig. 3) . 
Each vertex also has physical states. The number of the 
states depends on the states on the connected edges and 
they are labeled by a — 1, ...,Nijf. (see Fig. 2). Here we 
assume that 



N, 



ijk 



N 



jki- 



(9) 



Despite the similar look between a index and a index, 
the two indices are very different, a index labels the 
vertices while a index labels the state on a vertex. In 
this paper, we very often use a to label states on vertex 
a. 



The states on the edge are always bosonic. However, 
the states on the vertices may be fermionic. We introduce 
Sijfc(cc) to indicate whether a vertex state labeled by a is 
bosonic or fermionic: Sijk(a) — if the a-state is bosonic 
and Sijk{o) = 1 if the a-state is fermionic. Here the 
vertex connects to three edges i, j, and k (see Fig. 2). 
The function s,jfc(a) satisfies 



(10) 



Only a fermionic operator can map a bosonic state to a 
fermionic state and vice versa. Let N^ k {Nf- k ) be the 
number of bosonic (fermionic) state on the vertex that 
connect to three edges with states \i), \j), and \k). We 
also have 



N h l3k 



N 



f 



ijk 



N 



f 



(11) 



Each graph with a given a,/3, j, ... labeling (sec Fig. 
2) corresponds to a state and all such labeled graphs 
form an orthonormal basis. Our fixed-point state is a 
superposition of those basis states 



all conf. 




(12) 



In the string-net approach, we made a very strong as- 
sumption that the above graphic states on two graphs are 
the same if the two graphs have the same topology. How- 
ever, since different vertices and edges are really distinct, 
a generic graph state does not have such an topological 
invariance. To consider more general states, in this pa- 
per, we would like to weaken such a topological require- 
ment. We will consider vertex-labeled graphs (v-graphs) 
where each vertex is assigned an index a. Two v-graphs 
are said to be topologically the same if one graph can be 
continuously deformed into the other in such a way that 
vertex labeling of the two graphs matches. In this paper, 
we will consider the graph states that depend only on the 
topology of the v-graphs. Those states are more general 
than the graph states that depend only on the topol- 
ogy of the graphs without vertex labeling. Such a gener- 
alization is important in obtaining interesting fermionic 
fixed-point states on graphs. 



B. The structure of a fixed-point wave function 

Before describing the wave function renormalization, 
we examine the structure of entanglement of a fixed point 
wave function ^g x on a v-graph. First let us introduce 
the concept of support space with a fixed boundary state. 

We examine the wave function on a patch, for exam- 

' j k / i j k\ 

pie, ■ The fixed-point wave function f/;fi x ( j 

(only the relevant part of the graph is drawn) can 
be viewed as a function of a,f3,m: (pijkl,r( a i A m ) = 



6 



'\ J. k\ 



if we fix i,j,k-,l and the indices on the 



other part of the graph. (Here the indices on the other 
part of the graph is summarized by T.) As we vary the 
indices T on the other part of the graph (still keep i, j, k, 
and / fixed), the wave function of a, (3, m, </>ijfcz,r(a, P, m), 
may change. All those 4>ijki,r{ct, f3, m) form a linear space 
of dimension Dijki*. Dijki* is & n important concept that 
will appear later. We note that the two vertices a and j3 
and the edge m form a region surrounded by the edges 
i,j,k,l. So we will call the dimension- -Dyjw* space the 
support space Vuki* and D^ki* the support dimension 
for the state ?/>fix on the region surrounded by a fixed 
boundary state i, j, k, I. 

We note that in the fixed-point wave function 



by tpfi-x \§^n • Again, such a wave function can 
be viewed as a function of x, S, n: (ftijki r(Xi ^ n ) — 
ipRx if we fix i,j,k,l and the indices on other 

part of the graph. The support dimension of the 
state ipa x on the region surrounded by k, I 

is Dijki* ■ Again D ijk i* < N ijkl *, where N ijkl * = 
12n=oNkjn*Ni» n i is number of choices of x>&-> n - The 
saturation assumption implies that Nij k i* = D^u* ■ 

The two fixed-point wave functions ifjRx I j^v I and 



We note that in the fixed-point wave function The two fixed-point wave functions Vfix anc 

(u }\ V ' / 

"0fix „\K , the number of choices of a, p,m is /' 7 oc/\ 

\ 1 J ipRx \§^n are related via a local unitary transforma- 

Nijki* = J2 m =o N jim*N km i*. Thus the support dimen- V ' / 



I2m=o N 3im* N kml*- Thus the support dimen- 
sion Dijki- satisfies Dijki* < Nijki* ■ Here we will make 
an important assumption - the saturation assumption: 
The fixed-point wave function saturates the inequality: 



D^ki* = Nt 



ijkl* 



N 

^ Nji m *Nkml* 



(13) 



tion. Thus 
which implies 

N 



Dijki* Dijki* ) 



N 



l*ni- 



(15) 



(16) 



In general, we will make the similar saturation assump- 
tion for any tree graphs. We will see that the entangle- 
ment structure described by such a saturation assump- 
tion is invariant under the wave function renormalization. 
Similarly, we can define Dijk as the support dimension 



m=0 



71=0 



We note that the support space of -0fi x | jxP^ I and 



of the $fi x 



on a region bounded by links i,j,k. 



Since the region contains only a single vertex a, we have 
Dijk _• Nijk- The saturation assumption requires that 



D 



ijk 



N, 



ijk- 



(14) 



In fact, this is how N^ k is defined. 

We note that under the saturation assumption, the 
structure of the support dimensions for tree graphs is 
encoded in the Nij k tensor. Here Nijk plays a similar 
role as the pattern of zeros in a classification of fractional 
quantum Hall wave functions. 35 



C. The first type of wave function renormalization 

Our wave function renormalization scheme contains 
two types of renormalization. The first type of renor- 
malization does not change the degrees of freedom and 
corresponds to a local unitary transformation. It corre- 

> A * ' J y k 

sponds to locally deform the v-graph to . 

(The parts that are not drawn are the same.) The 
fixed-point wave function on the new v-graph is given 



V'fix N^n should have the same number of fermionic 
states. Thus eqn. (16) can be splitted as 

N 

rn — 
N 

= £(iV^iV£m + N[ jn . Nf.J, (17) 



N 



J2( N !irn*N( ml .+Nf im ,N b kml *) 
N 

Y^{Nt jn ,N(, m + Ni jn ,N^ m ). (18) 



We express the above unitary transformation in terms 



of the tensor F^™^, where a, b = ±, i,j,k,.. 



kln : x8,b ' 

0, N, and a = 1, Nijk, etc : 

N N kj „* N nil , 

bijki,T{a,P,m)~Y^ 12 F kin,xS,xt^Jki.r(xJ,n) 

n=0 x=l <5=1 

(19) 
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or graphically as 



fix 



E ij'm, a/3, o/3 
I 'kln,x5..x5 ^" fix 




(20) 



where the vertices carrying the states labeled by 
(a, (3, x,S) are labeled by (a, (3, %, 6) (see Fig. 2). Here 
~ means equal up to a constant phase factor. (Note that 
the total phase of the wave function is unphysical.) We 
will call such a wave function renormalization step an 
F-move. 

In general, we should write the F-tensor as F kln xS '~' s ~ 
which depends the vertex labels a, j3, x, 5. In this paper, 
we assume that the F-tensor only depends on the sign of 
P — a. It does not depend on how big is the difference 

I n I <n 7-1 • 7-,^ m j a /3> a P 

\p — a\. So the b -tensor can be written as b kln xS x j~. 
Here a/3 is a function of a and (3 describing the sign of 

[3 - aT 



af3 



if a > (3, 



+ if a < 13. 



(21) 



There is a subtlety in eqn. (20). Since some values of a, 
/3, ... indices correspond to fcrmionic states, the sign of 
wave function depends on how those fermionic states are 

ordered. In (20), the wave functions -0fix ^^^P 7 ^ I an d 

&n are obtained by assuming the fermionic 
states are order in a particular way: 



fix 




aj3rjir]2. 



= > A ■ 



W5mm-)x&n ,n ... (22) 



where sum over the indices on the vertices and on 
the edges and rji are indices on other vertices. Here 
\aj3...} otfi... is a graph state where the a-vertex is in 
the |a)-state, the /3-vertex is in the |/3)-state, etc . 
We note that \(3a...) p a ... is also a graph state where 
the a-vertex is in the |a)-state, the /3-vertex is in the 
|/3)-state. But if the |a)-state and the |/3)-state are 
fermionic (ie s jim *(a) = s fcm/ »(/3) = 1), \a(3...) ^p... and 
\/3a...) f3 a... will differ by a sign, since in \a/3.-) a p... the 



fermion on /3-vertex is created before the fermion on a- 
vertex is created, while in \0a...) a/3... the fermion on a- 
vertex is created before the fermion on /3-vertex is cre- 
ated. In general we have 

\a0...)«p... = (-) s -"*^ Sfc -*^|/3a...)^.... (23) 



We see that subscript a/3... in |a/3...) a /3... is important 
to properly describe such an order dependent sign. Sim- 
ilarly, we must add the superscript in the wave function 

as well, as in ij)^- 1 - 2 j > since the amplitude 

of the wave function depends both one the labeled graph 
! „ 'l k 

and the ordering of the vertices afif] r] .... Such 
a wave function has the following sign dependence: 



, aPv rj ... 
V>fi x 



= i_y jtm . (a)a fcmI . GSty£ s 2, V" 



(24) 



Thus eqn. (20) should be more properly written as 



fix 




E ijm,aP,ap x8.. 
fcZn,x<5, 'fix 



(25) 



where the superscripts a/3... and x£--- describing the 
order of fermionic states are added in the wave function. 

Since the sign of the wave function depends on the or- 
dering of fermionic states, the F-tensor may also depend 
on the ordering. In this paper, we choose a particular 
ordering of fermionic states to define the F-tensor as de- 
scribed by a/3... and x<5... in eqn. (25). In such a canon- 
ical ordering, we create a fermion on the /3-vertex before 
we create a fermion on the a-vertex. Similarly, we create 
a fermion on the 5-vertex before we create a fermion on 
the x -ver tex. 

In the following, we would like to introduce one Ma- 
jorana numbers 9 a on each vertex a to rewrite (25) in 
a form that will be valid for any ordering of fcrmionic 
states on vertices. The Majorana numbers satisfy 

B\ = 1, 9a0p = -OpOa for any a ^ /3, 

f s = 0a, {e & ...etf = o t ..9 & . (26) 

We introduce the following wave function with Majorana 



numbers: 



lix 



[e%* 



\4> 



iix 




(27) 

where the order of the Majorana numbers 9 a 9 «... is tied 
to the order a/3., in the superscript that describes the 
order of the fermionic states. We see that, by construc- 
tion, the sign of does not depend on the 
order of the fermionic states, and this is why the Ma- 
jorana wave function ^^m^ 7 ^^ does not carry the 

superscript a/3... 

We like to mention that (0a,9p) and (6 x ,9s) are 
treated as different Majorana number even when, for ex- 
ample, a and x take the same value. This is because 
a and \ label different vertices regardless if a and \ 
have the same value or not. So a more accurate notation 
should be 



T, 4>i]ki,v{x^,n) and <f>ijki,v(a, fi,m) span two N ijk i> di- 
mensional spaces. Thus we require that F^^j to be an 
Nijki* x Nijki* unitary matrix 

E ijm'.a'^', a.p . ijm,a/3, _ , . 

^fc;«,X<5,X^ \ t kln,x&, X^ > — °rnaP,m'a'P'i W 1 J 

nxS 

where &map,m'ot'P' — 1 when m — ml ', a = a', /3 = /3', 
and S m a.pm'a'P' = otherwise. This way, the F-move 
represents a fLU transformation. 

We can deform to "1^, and to 

We sec that 



7 / nxS 



fk I O 



(32) 



Eqn. (30) and eqn. (32) relate the same pair of graphs, 
and should be the same relations. Thus 



-r-ijm,a./3 .a j^kl* m* ,fta, — a 



(33) 




i j v k\ 



'i 



(28) 



where 9 a and 9 X are different even when a = x- -But in 
this paper, we will drop the ~ and hope that it will not 
cause any confusions. 

Let us introduce the F-tensor with Majorana numbers: 



F, 



ijm,ap,a 



= (\s^ m *(a)s kml *(l3)+s kjn ,(x)s n ii*(S) Tpkl*m*,0a-a ,oa\ 



(where we have used the saturation condition Nijki — 
Dijki and a = a/3 — — /3a.) 

Using the relation (31), we can rewrite eqn. (25) as 



mat.fi 



'' J. 



(35) 



.ijm,aP,aP _ flSjim , (a) flSfcm! , {P) a s nil , (S)^s kjn , ( X ) w iim,aP,a§_ 



T' 

■ r kln,x&, x$ ~ - P 

We can rewrite (25) as 



F, 



or as 



(29) 



lix 



E ijm,ap, o/3 



nx$ 



(30) 



Such an expression is valid for any ordering of the fermion 
states. ., 

For fixed i, j, k, and I, the matrix F^'—g with ma- 

trix elements (F fei > ^) n ^i = F kin, x s,x]r 1S a matrix of 
dimension N^i* (see (16)). Here we require that the 
mapping 4> ljk i t r{x,S,n) -> cf) ijkLr (a, /3,m) generated by 
the matrix -j- to be unitary. Since, as we change 



1 J. k \ 



\ '/ / ma/3 \ 1 / 

(36) 

where 



/ -r-ijm, a/3, a\f 



~ d X d S_ V p Vol {^kln, x S,b ) ■ ( S7 > 

Here we have used 

(a)^ m! . (/3)^ n<! , (x) )t 



/jSJsJn* (x)/jS n ii* (5) fl s fcmi , (P) a s jim * (a) 
- y x tf «5 61 /3 y <3; 



(38) 
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We can also express ^fi x ( j as 



fix 



' 1 1 %} 



maft 




using a relabeled eqn. (30). So we have 



(T- 

W kln, X S,b 



■pjkn,xS,b 
l*i*m* ,/3a, — a ' 



(39) 



(40) 



The above relation between ^ 



fix 



<i 



and 



^fix ( should be the same regardless if we choose 



X6 
E 




Therefore, we have 



■j~ijm,afi,a / -r-ijm' ,a' j3 f X ^ -r~ijm,a/3,a / -pijm' ,cx' 0' ,b\-\ 

J-ui„ „x , k- / "fe«n,x5,+ > — / ^kln.-vS - \ J 'kln,xS,- > 



kln,xS,+ 



(46) 



Since the total phase of wave function is unphysical, we 
can adjust the total phase of -^™^f b ' a for each choice of 
(a, b) — (±,±) independently to make ~ into =: 



/ -pijm,a.0 .a\f -r-jkn,x$,b 

\ kln,x<>;b ) •'l*i»m* ,/3a, — a' 



\ r kln,xS,b I ~ V ) 



However, the following phase change 



ijm,a.l3 ,a 



(41) 
(42) 
(43) 



does not affect eqn. (42) and is still not fixed. 

If we apply eqn. (42) twice, we reproduce exactly 
eqn. (34). Thus eqn. (34) is not independent and can 
be dropped. If we apply eqn. (42) four times, we get the 
identity transformation. Applying eqn. (42) once, twice, 



F, 



and three times allow us to express F k j n s _ , ± kin x $ + , 



F 



in terms of F, 



ijm,a0 ,4- _ 



fcin,x<5,+ ' <-cj.ij.jo ui -i kln,xS,- 

T?ijm,aP,+ _ (\s jim *(a)s kml *(j3)/-pjkn,xS,+ \* 
r kln,xS,+ ~ V ) \ r Z*i*m* ,0a,-) ' 

pijm,aP - _ r_y Jlm *{a)s kml ,{p) + s kjn *( x )s nil *{5)pkl*m*,[3a,+ 
kln,x&,+ ^ ' ij*n* ,6x>— ' 

rpijm.afi _ (_\s hjn * (x)s„a* (*) f pJ*inVx,+ \* (AA\ 
r kln,x&,- ~ \ > \ r jk'm,aj3,~) ' 

We can use transformation eqn. (30) and eqn. (36) to 
rearrange the labels of the vertices: 



\ U j n X S \ n 



ijm,aj3, a/3 , ijm f ,a' fl' , a'/3' . + 



nx<5;m'a'/3' 



T 

kln,xd, x$ 



(•^Un,xS, 



XS 



(45) 



Since a/3 and a'j3' can be different, the labeling of the 
vertices is rearranged. 



r kln,xS,+ ^Wn,xt5,+ J 

nx<5 



iSJ. 6 T?ijm,aP,a / -rpijm' ,ct' 
e 2_^ r kln,x&,- \ r kln,x&,- 
nxS 



(47) 



The above conditions on the F-tensor ensure the self con- 
sistency of rearranging the labeling of the vertices. 
From the graphic representation (25), We note that 



F, 



N jim * < 1 or N km t* < 1 or N kjn * < 1 or N nil * < 1, 

or Sji m * (a) + s km r* (/?) + Sfcin* (x) + (<5) = 1 mod 2. 

When Nji m * < 1 or Nkmi* < 1, the left-hand-side of 

eqn. (25) is always zero. Thus F^[ a f^ a = when 
AT, 



< 1 or iVfc m i, < 1. When N kjn * < 1 or iV m; . < 1, 
wave function on the right-hand-side of eqn. (25) is al- 
ways zero. So we can choose F$™£fi} a = when N k j n * < 

1 or N ni i' < 1- Also, F % k ^^^ a represents a pseudo- 
local bosonic operator which contains even numbers of 
fermionic operators. Therefore F k l j^ l ^f^ a is non-zero only 

when sjim* (a) + Skjn* (x) + smi* (S) = mod 

2. 

The F-move (30) can be viewed as a relation between 
wave functions on different v-graphs that are only differ 
by a local transformation. Since we can locally trans- 
form one v-graph to another v-graph through different 
paths, the F-move (30) must satisfy certain self consis- 

i j k I 



tent conditions. For example the v-graph m \0v c<ni 

p 




be transformed to 




through two different paths; 



one contains two steps of local transformations and an- 
other contains three steps of local transformations as de- 
scribed by eqn. (30). The two paths lead to the following 
relations between the wave functions: 
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fi J k l \ 



lix 




E ijm,ap, up 
knt,r}tp,r]Lp h x 




E ijm,a@, a/3 itn.ipx, 
knt,r]tp, r]ip lps,n^,Kj ^ x 




E ijm.a/3, aft itn.ipX: <PX j~jkt,7]K, V K ^j 
knt,r)<p,rjLp lps,K~y, K~y lsq,5(fi,5(fi ^ x 




(49) 



lix 



A j k l\ 

V p ) 



E mkn,t3 X ,l3x r 

q5e 




E 

q8e;s<f>j 



mkn,0x, Px T-ijm,ae, ae, T , 



lpq,5e, r5e 



qps.(pj, ^>7 



V 



/' 



(50) 



/ 



The consistence of the above two relations leads a condi- 
tion on the f-tensor. 

To obtain such a condition, let us fix i, j, k, I, p, 



and view V'fix 




as a function of a, j3, x, m, n: 



As we vary indices 



on other part of graph, we obtain different wave func- 
tions (j)(a, /?, x, 77i, n) which form a dimension Dijkip* 
space. In other words, Dijkip* is the support dimension 
of the state V'fix on the region a, /?,%, m, n with bound- 
ary state i,j,k,l,p (see the discussion in section VB). 
Since the number of choices of a, j3, x, m, n is Nijkip* = 
I2 m ,n N 3 im* N kmn * N inp * , we have D ijMp * < N ijk i p *. 
Here we require a similar saturation condition as in (13): 



Nijkip* Dijkip* 



(51) 




Similarly, the number of choices of S, (f>, 7, q, s in 



is also Nijkip* ■ Here we again assume 



Dijkip* = Nijkip- , where Dijkip* is the support dimension 

Of t/>; 



on the region bounded by i, j, fc, 



So the two relations (50) and (49) can be viewed as two 
relations between a pair vectors in the two Dijkip* dimen- 
sional vector spaces. As we vary indices on other part of 
graph (still keep k, l,p fixed), each vector in the pair 
can span the full Dijkip* dimensional vector space. So the 




I 

validity of the two relations (50) and (49) implies that 

N hit* N Un , N lt3 * 



t >] — 1 if — 1 K — 1 

N * 



ijm,a/3, a/3 itn,ipx, V?X - 1 -j^: r ) K --: 
knt,r}(p, r)ip Ips^n-y, hey lsq,5<fi, S(p 



qps,cj>-y, (j>-y ' 



^~lpq,5e, 5e **~ i 



(52) 




which is the fermionic generalization of the famous pen- 
tagon identity. The above expression actually contain 
many different pentagon identities, one for each labeling 

i j k I i j k I i j k I 

scheme of the vertices in m \z/ , nt \e/i , , 

p p p 

i j k l i j k I 

X%v 

, and \A ■ Since the vertex labeling can 
p p 
be rearranged through eqn. (45), so we only need to con- 
sider one pentagon identity from a particular labeling 
scheme of the vertices: a="/=Lp<(3=e=(f)=ri< 
X = 6 = k. We obtain 

N kjt* N Hn , N Hs , 

EV * V * \ ^ j~ijm,a{3.+ -pitn.^px^ -pjkt,riK,, J r 

/ j / j / j knt,rjip,— lps,K*y,— lsq.&4>>~ 
t ri—1 ip—1 k—1 

N * 

— Z-^i ,r lpq,5e,- ^qps,^- ■ V JO ) 

6=1 

We can use the transformation (43) to change ~ to = in 
the above equation and remove the Majorana numbers 
to rewrite the above as 

N kjt * N Un . N lt3 * 

EV 1 \^ pijm,aP,+ pitn,tpx,+ pjkt,t)K,+ 

/ j / j / j knt,r\ip,— lps,K~(,— lsq,5ij),— 
t rj—1 ip—1 k—1 

N 

= (_ynm* (a)s(fcg* (<5) \^ pmkn,0x,+ pijm,ae,+ ,r.\ 
V / / j lpq,8e,— qps,07,— ' \ ) 
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The above fermionic pentagon identity (54) is a set 
of nonlinear equations satisfied by the rank-10 tensor 
The above consistency relations (44), (47), 



^ijm,a8,-\- 



and (54) are equivalent to the requirement that the local 
unitary transformations described by eqn. (30) on dif- 
ferent paths all commute with each other up to a total 
phase factor. 



D. The second type of wave function 
renormalization 

The second type of wave function renormalization does 
change the degrees of freedom and corresponds to a gen- 
eralized local unitary transformation. One way to im- 
plement the second type renormalization is to reduce 

'» so that we still have a trivalent graph. 

This requires that the support dimension Dai* of the 



fixed-point wave function ip^ 



is given by 



5u> ■ 



This implies that 

V'fix 




(55) 



(56) 



The second type renormalization can now be written as 
(since D^* = 1) 

where the ordering of the vertices is described by a(3rj.... 
We will call such a wave function renormalization step a 
O-move. Here 0{ k ^ a ^ satisfies 

N kii , N rjk * 

^ ^ ^ Oi% a0 (Oi^T = 1 (58) 



k,j a=l 8=1 



and 



Oih.afi = Q) if N , km , t < 1; Qr N , fjk < ^ (5g) 

or Sik*j*(a) + 8i*jk(J3) = 1 mod 2. 

The condition (58) ensures that the two wave functions 
on the two sides of eqn. (57) have the same normalization. 
We note that the number of choices for the four indices 
(j, fc, a, 0) in Oj k ^f must be equal or greater than 1: 



D: 



jk 



*Ni* jk > 1. 



(60) 



In fact, we should have a stronger condition: the number 
of choices for the four indices (j, k, a, (3) that correspond 
to bosonic states must be equal or greater than 1 



Di 



jk 



N 



f 



nLu > i. 



ik*j* i* jk 



(61) 



The wave functions in eqn. (57) is defined with re- 
spect to the ordering of the fermionic states described by 
aj3rj.... Let us introduce 



*fix (4. 




and 



-,jk,a8 _ n s ik* j' ( a ) n s i* jk(P) n jk,a/3 
i, a/3 — ° 8 U i,a8 • 



We can rewrite (57) as 




jk,ap v f 



(62) 
(63) 

(64) 



which is valid for any ordering of the fermionic states. 
Notice that 



fix 




^<^ Q * fix (-U). (65) 



Since eqn. (64) and eqn. (65) are really the same expres- 
sion, we have 



(66) 



Since the total phase of the wave function is unphysical, 
we can adjust the total phases of O^Y 1 ^ and oj^"^ in- 
dependently. We can use this freedom to make eqn. (66) 
into Oj a p = eV,-a or 

Qjk,a/3 _ /_y lk * 3 ,(a)s jki ,(B)Qk*j\8a^ 



O 



We can find another relation between O^f 1 ^ and 

jk,a8 



E ij*k,al3,af3_ 



E ij"k,a/3, afl —j*jm,vn, v/j. 
jim,vfi, u^i ^"i*i*n*,A7, A7 fix 




~> ■ _ ii/a J > 

(68) 

which leads to 

/njk,u8 ^ ST^ -pij*k,a8,a -pj*jm,ufj,,c snjn,y\ / fi q\ 
i,a — / j • r jim,vfi,c i*i*n*,A7,-6 i,b ' 

After removing the Majorana numbers, we find 

Of a ' a " ^ (70) 

\ ) jim.y^.c i*i*n* ,X*y, — b i,b 

rrifiw^wyX 
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Choosing a — — and b = c = + in the above, we find 
another relation between Oj 1 ^ 01 ^ and O^j ^ . Combining 
it with eqn. (67), we find 

E(_yi„* } * h)s Jni * (A) pij'k,af!,- rf*jm,^,+ 
\ / jim,u[i,+ i*i*n* ,A7,— i. 



n,7A 
+ • 



tained by noticing 



iix 




E. The fixed-point wave functions from the 
fixed-point gfLU transformations 

In the last section, we discussed the conditions that 
a fixed-point gfLU transformation {F^^ff^ 1 , 0{ k ^ a ^) 
must satisfy. After finding a fixed-point gfLU transfor- 
mation (Flj™^ 3 ^, Of k ^ a ^), in this section, we are going 
to discuss how to calculate the corresponding fixed-point 
wave function ipfi K from the solved fixed-point gfLU trans- 

formation (i^^<\0|^). 

First we note that, using the two types of the wave 
function renormalization introduced above, we can re- 



^2 ^Mnx* x*^ &x ( J dUCC graph t0 Cj)- So ' once we know ^ fix (G^) ' 



More conditions on Fj,j™^f^ a and 0{ k ^ a ^ can be ob- 



"XA 



E mjp,ri/3, riP 
m,aA, a A j 

"XA;97<5 




we can reconstruct the full fixed-point wave function V>fix 
on any connected graph. 
Let us assume that 



iix 



= A 1 



(78) 



which implies that 



Qjm,ar, s ^ 



(73) 



/ 



■^kln^-.x^ nlk,-f5,~fS U fe*,7X Wfix 



i V 

\5/ 



nxA;7<5 

We find that 

/ j kln,x^> xA nlk,-y5,-y8 &* ; 7X <5 
"XA7 

for all k,i, I satisfying Nku- > 0. 



(74) 



If we choose a particular labeling of the vertices 

a = S<rj = X = ^/<f3 = x, (75) 
we reduce the above expression to 

r^, ST^ -p™3P,vfl,+ 7ry* m . Q A, + ^,jn*,7X/)S fcil » (S) 
— 2-^1 W",xA,- %r nlk,iS,- ^k*,+ U S 
«XA7 

for all k, i, I satisfying Nku* > 0, (76) 



e ie °iOi™' ar, 5 ip 5f,s = {-) Sk ^' ix)[Sk ^' ix)+s ^ iv)] x 

T?™3P,V@,+ Tpij*m,a\,+ n jn* ,fx 
r kln,x\~ r nlk,fS.- W k*,+ 

for all k, i, I satisfying Nku» > 0. (77) 



Note that "I^x ^(^T)^ contains no Majorana numbers 
and thus is equal to ^g x ^(^7)^ • Here A 1 satisfy 

J2A l (Ay = l, (79) 



which is simply the normalization condition of the wave 
function. 

To find conditions on A 1 , we note that 



This allows us to show 



or 



o( a k *' a0 A l * = e i9 *0$*f' af> Af. 



(80) 



(81) 



(82) 



Also, for any given k,a that satisfy Nkji > 0, the 

must be non-zero for some 



wave function i/>fi x [ "H*- r ] 

V ^ J 



r, where T represents indices on other part of the graph. 
Then after some F-moves and O-moves, we can reduce 

r to a (^}' 3 - So, for any given i,j,k,a that sat- 
isfy Nkji > 0, ipsx ( a (^)^j is non-zero for some j3. 
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Since such a statement is true for any choices of ba- 
sis on the vertex a, we find that for any given i,j,k,a 
that satisfy N k j$ > and for any non-zero vector v a , 

^2 a v 'aip 'fix ^™(^)^ i s non-zero for some j3. This means 
that the matrix Mf.ji <a given by 



{Mkji, a )ap = Vfe ( a (^) P ] = OQ'^A 1 * (83) 



is invertiblc: 

det(M fciil0 ) £ 0. (84) 

The above also implies that 

N kji = N i . j . k .. (85) 

Such a condition ensures that eqn. (60) is always valid. 
So we will drop eqn. (60). Due to the mod 2 conservation 
of fermion numbers, eqn. (59) also allows us to show that 



N, 



f 



(86) 



kji k* ' 

If we further assume that the wave function is well 
defined on a sphere, we can have additional conditions. 

One of them is ipf^ ^(^T)^ = V'fix (C^)\ which leads 
to 



A 1 = A l 



(87) 



The second one is A 1 ^ 0. To show this, let us first 
consider the fixed-point wave function where the index on 



a link is i: "0fix(*, T) = ipfix ^ T ^L- j , where T are indices 
on other part of graph. We note that the v-graph r _ 

i 

can be deformed into the v-graph r on a sphere. 
Thus V'fix ^ r — '— ^ = "0fix ^ ^ • Using the F-moves 
and the O-moves, we can reduce 



to Q: 



lix 



r i 



lix 



lix 



We see that 



(89) 



for all i. Otherwise, any wave function with i-link will 
be zero. 



* fix | a ^P), which implies OQ'^A 1 ' 

k : a 



The third condition comes from 



Ol i ' a/3 A fc * or 



Qj'k" ,Q/3 j^' _ qC j* ,ap ^k* 



(90) 



In fact the above relation and eqn. (79) allow us to show 
9 a A = in eqn. (82). When 6\ — 0, above relation can 
be derived from eqn. (79) and eqn. (82) and thus can be 
dropped. The conditions (79, 82, 84, 87, 89) allow us to 
determine A 1 . 



F. Summary of fixed-point gfLU transformations 

To summarize, the conditions (16, 18, 61, 85, 86, 31, 
42, 47, 48, 54, 58, 59, 67, 71, 77, 79, 82, 84, 87, 89) 
form a set of non-linear equations whose variables are 
N vk, N/ jk , Fl 1 ™^/, Of^" 13 , A\ and the seven phases 

(e ie ^\ a =±,b=±, e i0oi ,e ie ° 2 ). Let us collect those condi- 
tions and list them below 

N N 



m=0 
N 



n=0 



J2( N ^Nl ml ,+Nj m ,N b kml ,) 



m=0 
N 



• N kji = Nw , • = N^, k , . (91) 



E ijm ,a P ,aP . ijm,af),ul3. t _ . 
*"kln,x&,X§_ \ kln, X S,xl > ~ °maP,m'a'p' , 

/jpijm,aP,a\* _ (_\Siim* (°O a kml* (0) pjkn,xS,b 

\ kln,xS,b I \ J l*i*m* ,Pa,~a' 

ET7<ijm,ap,a / rpijm' ,q' P' 
^kln, X S,+ \*kln,x$,+ I 



1X& 



A9 a p T?ijm,aP,a / rpijm' ,a p' ,b\* 

e Z_^ r kln,x&- y r kln,xS,- I > 
n\5 



• F, 



,ijm : a{3.a 



kln,x8,b = wnCn N jim* < 1 Or N km l* < 1 Or 

N kjn , < 1 or N n a* < l,or 

(a) + s km i*(/3) + Skjn* (x) + S n u. (S) = odd, 

Nkjt* N tin * N lts , 

Ey \^ pijm,aP,+ T?itn,tpx,+ pjkt,r)K,+ 

/ j / j / j knt.rjip,— /ps,K7,— lsq,S(p,— 



t Tj—1 <p=\ K— 1 



N„ 



r lpq,Se- r qps,4>i,- " \^ A > 
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j\[ k . . * Nj * j & * 

• EE E oj^(oi^y = i, 

k,j a=l 0=1 

• = 0, if N ik .j. < 1, or JVi. ifc < 1, 
or 8ik*j*(a) + Si*jk(/3) = odd, 

• e ie ° 2 O k y+ ,l3a = (-) s ^*j*( q ) s j<"*(/3) x 

E(_yin*i* (7)Sj„i* (A) pij*k,ap,- j^j*jm,u^,+ Qjn^X 

m/if;n7A 

• e lfloi Q££ ,a %,<S /w = E (-) Sfc -* (x)[Sfc -* (x)+s ^* ( " )] x 

^fcin.xA,- r nlk,yS- w k' ,+ 



for all fc, z, Z satisfying iV^y* > 0. 



det(M feii!Q/ 3) ee det 



(93) 




7^ 0. (94) 



Finding N ijk) Nf jk> F%™ffi, , and A* that sat- 

isfy such a set of non-linear equations corresponds to 
finding a fixed-point gfLU transformation that has a 
non-trivial fixed-point wave function. So the solutions 
{ N i3^ N lj^ F kln^\,b^ Ol k a a0 ,A l ) give us a characteriza- 
tion of fermionic topological orders (and bosonic topo- 
logical orders as a special case where N^ k = 0). 



VI. GENERALIZED TURAEV-VIRO 
INVARIANCE AND ITS GRAPHIC 
REPRESENTATION 

As we know, the string-net model has a very deep rela- 
tionship with the Turaev-Viro invariance 36 for closed tri- 
angulated 3-manifolds. 37,38 In this section, we will show 
how such a states sum invariance is generalized to its 
fermionic version. 



A. The amplitude of tetrahedra graph and 
fermionic symmetric Q symbol 

Let us define the amplitude of the tetrahedra graph 
in Fig. 4 as G^in'-xf ' We call such a Majorana number 
valued tetrahedra amplitude as the fermionic symmetric 
Q symbol if it is well defined on a sphere and has the 
twelve tetrahedral symmetry. Notice here we assume the 
Q symbol is independent of the labeling scheme a/3x<5. 
Such an assumption is a consequence of the tetrahedral 
symmetry and will be proved later. 




_ V~> q~ijm,afl,aP j 
K ~ 2-1 kln,x'S' 




8' n X 



X'S' 



FIG. 4. We define the Q symbol as the amplitude of the 
tetrahedra graph, which can be expressed in terms of J-, O 
and A\ 



In terms of Majorana number, we can express the am- 
plitude of the tetrahedra graph as: 



y kln,x$ 



E i T m,a0,al H * >S > S j* fc* , x > \ 

• r kln,x'6',x'S' n,6>S w n*,x'x 
X'S' 

Va Vp &x X 

E ij m , a f),aJl u* iS > S f ' k \x 'X An* 
r kln,x'S',X'$' U n,V6 U n*,x^x A 

X'S' 

U a U P U S U X ^kln,xS 



(95) 



with 
G 



ST p^™."/ 3 ^ n li*,S'S n fk*, X 'x An'(Qf{\ 
2^ t kln. X 'S\xrSy J n,S!S_ U n*,xbc A ^ b > 
X'S' 

Notice the above expression of Q makes it only depend 
on the labeling scheme of a/3 but not x'^'t x'x\ due 
to the relabeling equations on T and O. For example, no 
matter what value x8 takes, we can always chose x'x' 
■ . <ys = + and freely chose the value x'<5'. 

If the wavefunction is well defined on a sphere, then the 
Q symbol must satisfy the following twelve tetrahedral 
symmetry as shown in Fig. 5 (a),(b) and (c): 



^ kln,x8 



kl* m* ,(3a,/3a 



m ij ,a X ,ax /^-n 



The first equality is nothing but the consequence of 



T, 



kl*m* ,aj3,a{3 



kin,xs, x s ~ J ij*n»,6 x ,Sx , however, the second equality 



is a new condition on J- and O. 



The relation Q 



' kln,xS 



kln,xS 



,m 13 ,ax,ax 
^nklfiS 



implies that 



is independent on the labeling scheme a/3. 
This is because we can choose a/3 and ax equal to ± 
independently. Thus we can finally get rid of the label- 
ing a/3and rewrite the Q symbol as: 



-,ijm,a/3 _ fl «jim* («) flSfcmi* (0) (<5) Q s kjn» (x) Qijm,al3 



*kln,xS 



'0 



(98) 
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The tetrahedral relation can also be reexpressed as: 

s-*ijm,afl pkV 'm* .J3a s-tm*ij* 1 ax (QQ\ 

y kln, X S ~ y ij*n*,S X ~ y nkl,0S • V yy J 

The amplitudes of the corresponding mirror reflection 
of the above tetrahedron graphs (see Fig. 5 (e), (f), (g)) 
take different values but also have the twelve tetrahedral 
symmetries. 

y l*k*n',8x ~ y ji*n,xS ~ y n'k'l- ,0f 

B. Graphic representation for Q symbol 

1. Graphic representation in the bosonic case 

To explain the graphic representation of the Majorana 
number valued Q symbol, let us first consider a sim- 
ple case that all the vertex indices a/?x^ are bosonic 
(s jim *(a) = s km l*(0) = s„u*(S) = s k]n *(x) = 0). In 
this case, the Q symbol is equivalent to the G symbol. 

We can use a tetrahedra to represent the G symbol, 
the ten indices of the G symbol can be put on the six 
oriented links and four faces. The detail rule is defined 
as following: 

(a) We choose an arbitrary surface of the tetrahedra, 
put the three upper indices ijm of the G symbol on the 
edge of this surface and the lower indices kin on the rest 
three links of the tetrahedra. Notice for each upper index, 
the lower index should always be put on the opposite link 
on the tetrahedra accordingly. 

(b) We introduce a self-consistent rule to define the 
orientations of the links. We may put ijm onto the edge 
of the surface ABC with clockwise ordering and the ori- 
entations of i, j links are clockwise while the orientation 
of to link is anticlockwise(looking from the top vertex D 



FIG. 7. (Color online) The tetrahedra graphic representation 
of the G symbol is actually a due representation of the original 
tetrahedra graph. 

of the tetrahedra). The remaining three links all point 
to the top vertex D. Because three link indices in the G 
symbol uniquely determine a vertex index(e.g., the ijm 
link indices determine the a vertex labeled as a), it al- 
lows us to put the vertex indices of the G symbol onto 
the faces of the tetrahedra accordingly. Here we use the 
convention that the two upper indices a/3 correspond to 
faces ijm and klm, and the two lower indices x<5 corre- 
spond to faces kjn and iln. 

(c) If the orientation of a link i is opposite to the con- 
vention we choose in (b), then we simply replace it with 
its dual i*. 

Essentially, the above graphic representation is noth- 
ing but the dual representation of the original tetrahedra 
graph, as shown in Fig. 7. The twelve tetrahedra symme- 
tries of the G symbol are also automatically implemented 
in such a dual representation. 



2. Graphic representation for the full Q symbol 

Since the Majorana numbers only appear for those ver- 
tex indices, it is straightforward to generalize the above 
graphic representation into the fermionic case. We only 
need to put Majorana numbers 9g_ onto each surface of 
the tetrahedra to represent those fermionic states and 
have: 

r ijm,ap _ a s Jim*(o')QS km i*(0) a s na t(S) a s k j n *(x) r ,ijm,a0 
^kln, X S — U °L °p_ °5_ °X ^kln.xS • 

(101) 



16 



Notice here we order the Majorana number in the same 
convention as in the T symbol. 



C. The 2 3/3 2 move for the Q symbol 

1. The orthogonality condition 

To derive the 2 — > 3/3 — > 2 move for the Q symbol 
from the self-consistent equations of J- and 0, we need 
to assume the vertex indices in O satisfy the following 
orthogonality condition: 



O 



jk,a/3 



ik-j*{a-)a s i*ik{P) n jk c /ino-, 
I V/3 U i,aj3°a/3 (lU^J 




a 



p =0^ aP A' 

i ,af3 



Jy 




i X 




(b) 






This condition is automatically satisfied in the (bosonic) 
string-net theory and we believe it is also true for a large 
class of fermionic topologically ordered states. 

Under such an orthogonality condition, we can simplify 
the expression of G as: 



G 



kln,x$ 



A" 



(103) 



Because GkinxS 1S independent on the labeling schemes 
a(3, x'<5', 5^5, x'x> we can show F and O tensors are actu- 
ally also independent on those labeling schemes. Hence, 
we can omit the labeling indices for the Majorana number 
valued unitary transformation T and O: 



■ijm,a/3 



kln,x8 



'( a )nSkml* 
9 P 



n jk,a0 _ n s ik*j* ( a ) n s i*jk(0) n jk r 
W i,afi — °°L °/3 w i °<*P 



" (x) pijm,a/3 
kln,xS 

(104) 



Finally, we can omit all the labeling indices and sim- 
plify the expression of G as: 



G 



kln.xS 



pijrn,a/3 n li* n j*k* .n* 
t kln, X & U ™ °n> A 



(105) 



FIG. 8. The graphic representation of the amplitude of O 
graph and its inverse. 



See in Fig. 8 (b), using the right hand rule, the graphic 
representation indicates the cyclic symmetry as well as 
the 180 degree rotation symmetry: 



fpkji,al3 _ ip.fk'i* .J3a 
V a/3 — V j3a ' 



(108) 



which is automatically satisfied from the condition 

Similarly, we can use the left hand rule to define the 
inverse of ©(for those nonzero component) as: 



(e^y 1 = 9 s i im eT k ^ (a) (<;;;. r) l s ijk , 

(109) 



where dijk = 1 for N^k > and Sijk = for Nijk = 0. 

If O symbol satisfy the orthogonality condition, we can 
further simplify the above expression for graph as: 



2. Graphic representation for the & symbol 

To derive the 2 — > 3/3 — > 2 move for the Q symbol, we 
also would like to introduce the graphic representation 
for the 6 symbol, which is defined as the amplitude of 
the 6 graph. 

In the dual representation, the graph can be repre- 
sented as a triangle, see in Fig. 8. 



tpkji,aP _ n jk,afi A i 
V a/3 - U i*,a/3 Ji 



1 9L °0 U i',ap A ■ 



(106) 



We put the kji indices on the links of the triangle and the 
indices a/3 on the front surface and opposite surface of 
the triangle, then we use the right hand rule to determine 
the orientations of the links and the ordering of the a/3 
indices. If the 9 graph is well defined on a sphere, we 
can derive the cyclic symmetry: 



fckji,af3 _ r^ik^aP _ ,ry'i/c,a/3 
V a/3 — V aP — V a/3 



(107) 



(110) 



where the amplitude of graph is also independent on 
the labeling scheme. The inverse of (B can also be simpli- 
fied accordingly: 

= ^m^^yW (o?*a*) _1 6 a0 S ijk . 

(Ill) 

In this case, the cyclic condition (107) on symbol 
will further imply: 



OilfA* = 0%A j = Q%A k . 



The 2 -> 3/3 -> 2 move 



(112) 



Let us first consider a simple case that all the ver- 
tex indices af3\5 of Q symbol are bosonic(sji m . (a) = 
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Skmi*(f3) = s nil *(5) = Skjn»(x) = 0)- In this case, the Q 
symbol is equivalent to the G symbol, and the 8 symbol 
also has the simple form: 



From the pentagon equations of F: 



(113) 



Epmkn,0x pijm,ae 
lpq,Se qps,<p-/ 



= EE F < 

t 7],(p,K 



ijm 7 Oi/3 Tpitn,(px rpjUt^K 
knt,r)tp lps,K-y lsq,5<f) 

(114) 



we have: 



E 



U lpq,Se 



^iiym,ae 



o p q m * o L q : 1 * a*- o?*oil g *A s - 



EE 

t T},ip,K 



knt,r)tp 



siitn,ipx 



G 



lsq,S<f) 



or o t 



IT a^ ol 



Op"A s ' O s q 3 O k q l l 'Ai* 



Using the cyclic relation (112), we further have: 



E lpq,Se gps,07 
O p q m *Ai* 

lpq,8e qps,<pj 



fiijm,a/3 fiitn,ipX/~ijkt,r}K 
knt,rjip Ips.K^f lsq,8(p 

, ^OrOp'At'OP'A** 



2. 



O v q m Ai' 



EE 



A l G 



G 



knt,rj(f IpSjKj Isq 
m* A f Q rk' At p* A s ' 

(116) 



(115) 

In the last equation, we use the condition A 1 = A 1 . 

In terms of 8 symbol, we can further simplify the above 
equation as: 



E^<mkn,l3xr-fijm,ae' I r\m* pq" .ee'\ 
U lpq,8e U qps,^ { U ) 

ee 1 

= E E At ^;^GtT^ x <'f/ (e^'^y 1 (e**^)" 1 (<-) r ' -"') ' . (117) 



The generalization to the fermionic 2 — > 3/3 — > 2 move is straightforward and we only need to replace the G and 8 
symbol with their corresponding fermionic version: 



pq ,ee 



EQmkn,fixQijrn,ae' ( (D m *J 
^*lpq,5e ^*qps,<t>~( \ 

ee' 

E E ^^^C;i x C;i K ' (e^*^) -1 (e**^)" 1 (e^'^'y 1 . (ns) 



Notice for the fermionic CP 1 symbol, we need to apply 
the right hand rule to determine the ordering of the Ma- 
jorana numbers. Actually this is crucial to derive the 
correct fermionic pentagon equation for T: 



E-pmkn,Px -pijin.ae X ^ X ^ -pijm,a(3 -pitn,cpx -pjkt^K 
lpq,Se gps,07 / * / * knt,r]<p lps,n~f lsq,S(p ' 
e t rj,ip,K 



In terms of J 7 , the ordering is hxed because they do not 
commute with each other. However, in terms of Q symbol 
in eqn. (118), we don't need to specify the ordering since 
all of them carry different Majorana numbers, but we do 
need to put all the Q symbol in front of 8 symbol. 

From the graphic representation of the Q symbol, it is 
straightforward to represent the above generalized pen- 
tagon identity of the Q symbol as a two tetrahedron to 
three tetrahedron move(and its inverse), see in Fig. 9. 



D. The 2 ^> 0/0 -> 2 move for the Q symbol 

1. The associativity condition of A 1 

To derive the 2 — > 0/0 — > 2 move for the Q symbol, we 
need to further impose another associativity condition 
for A 1 : 

^25 k , l3 A l A j (-l) Sk '^^ = DA k , (119) 

ija 

where D is a constant. Notice the factor (— l) s fc*ij( a ) 
is a new feature uniquely arises from our fermionic gen- 
eralization. In the bosonic case, this equation can be 
simplified as: 

Y.foiiA'A* =DA k , (120) 

ija 

The loop value A 1 is proportional to the quantum di- 
mension, which should be a positive real number for 
most unitary theories. However, A 1 can be negative in 
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m. 




n 



FIG. 9. Graphic representation of the pentagon identity for 
Q symbol, Each tetrahedra represents a Q symbol. When 
two tetrahedral share a face, we insert the G _1 symbol as a 
metric and trace over both the front face and opposite face 
indices for that face. In this figure we omit the face indices 
for simplification. 



£A n A m A j 



nmj 




(-l^'D© 1 * 11 



FIG. 10. Graphic representation of the 2 — > 0/0 — > 2 move 
for Q symbol. In this figure we also omit the face indices for 
simplification. 



our fermionic generalization and we believe the factor 
(—iy><*ij( a ) play an important role to maintain the uni- 
tary properties of the theory. Such a conjuncture of the 
associativity relation on A 1 is consistent with other condi- 
tions but the underlying mathematical structure is still 
unclear, which will be studied in details in our future 
publications. 



In terms of G, we have: 
A 1 ' 

~ O l j" A k ' O 1 ™* A n ' 3 ,*, k * A n 
Use the definition of symbol, we further have: 

Si 



(122) 



J W Qtnji* ,r)otQi*lk* ,5/3 



"■XX 



(123) 



Combining with the associativity equation for A 1 , we 
can derive the 2 — > 0/0 — > 2 move for G symbol (see Fig. 
10): 



D Qi*lk*,5/3 _ y-] 

nm 3XX' A A' 77a 



An Am Aj^mji,riPr ii 3* m ' aX ' v 
Lr kln t xX Lr nlk,x*6 



r\m*ln* ,AA' 



Qmji ,77a 



(124) 



Similar to the 2 — > 3/3 — > 2 move, we can generalize 
the above expression to the fermionic case by replacing 
Q and with the corresponding fermionic version. How- 
ever, due to the anti-commutative relation between those 
Majorana number, a simple replacing does not work. Let 
us start from Eq. (74) and treat the sign carefully. Again, 
under the orthogonality condition of O, we have: 



Ej-mjp,r]P -pi]* m,a\ snjn 
• r kln,x>< • r nlk,x& k ' 

nxA 



(125) 



In terms of G, we have: 

or^ 5 u , k 8 ip 5 av 8 ps e s ^ (Q) ^-* ^ef*' {s) 



AnQmjp,r)f3 Qij* m,a\ 

, Of A k * Ol m * A n * Oll k ' A r 

nxA k n n 



E 



(126) 



2. The 2 -> 0/0 -»■ 2 move 



Recall the definition of (D, we can derive: 



To derive the 2 — > 0/0 — > 2 move, let us again start 
from the simple bosonic case. In this case, under the or- 
thogonality condition and tetrahedral symmetry condi- 
tion, we can simplify the J 7 , O relation Eq. (74) equation 
as: 



( l) s ""'-*( a )( l) s »*ifc* (P)lQ m 3i* ,r/a,gi*ifc* ,8/3 _ 

E An GZi p ;^G^f ) "" , '" xx .) ' (e m ' in *- xx ') 

"X-^X'^' 

(127) 

Notice the sign factors ( — l) s ^ji* (") anc [ (_]_)*<*i*»OT 
arise from the ordering of Majorana number. 
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Combining with the fermionic associativity relation of 
A 1 Eq. (119), we finally have: 



E 

-1 



An ™ Aj pmjitfP Mj* m,a\' 
71/1/1 y kln, X \ y nlk, X 'S X 



(®k*j*n,xx'^ ^Qm*ln*,xy\ ffo m ji*,r, a \ ^g) 

Same as in the 2—^3/3—^2 move, the ordering of Q 
symbol is not important as long as they appear in front 
of B _1 symbols. However, in terms of J-, O, the ordering 
is important and it is automatically taken care of by the 
ordering of Majorana number in (B _1 symbol, which is 
determined by the left hand rule. 



E. Topologically invariant partition function 

Based on the 3-^2/2^3 move and 2^0/0^2 
move, we can prove the following partition function is 
independent of the triangulation of closed triangulated 
3-manifolds. 39-41 



J top 



]JN V E E 
a/37<5--- ijklmn-- 



n w 

link tetrahedron face 



(129) 



where N v is the number of vertex. 

From the 2—^3/3—^2 move, it is obvious the above 
partition function is a invariance for those triangulations 
with the same number of vertices. To prove the invari- 
ance still holds for those triangulations with different 
number of vertices, we only need to show the correct- 
ness of 4 — > 1/1 — > 4 move based on the 2 — > 3/3 — > 2 
move and 2 — > 0/0 — > 2 move. The proof is very similar 
as the bosonic case, however, due to the fermionic na- 
ture, we need to take care of the sign. See in Fig 11, we 
want to prove: 

^ AnAmAjAt 

w' xx' 'fifi' KK> ' w' n f ri jt 

=^kln,x^ ^qst,K<fi ^tspjCp'o 1 =*qtp,pip 
X ^0*^* n *' K/t '^ 1 1 ^rn* st* ,<f>4>'\ 

X 1 ^© fc * 3* n ,xx'^ 1 fafirn* In* ,\\'\ 



DG, 



ikl,/3j 
qsp.pa 



(130) 



By applying the 3 — > 2 move, the left side of the above 
equation becomes 



\ " 



\ " 



-,mji,np pij* m,n' 'A' f->ikl,/3"f 



An Am Aj r mji,np r ij m,n a r% 
/ j / j A 71 71 y kln, X \ y nlk, X 'S y qsp,pa 

VV'XX 1 AA'55' nmj 



(131) 



£A n A m A j A' 

nmjt 




= D 



FIG. 11. Graphic representation of the 4 — > 1/1 — > 4 move 
for Q symbol. In this figure we also omit the face indices for 
simplification. 



Finally, we can apply the 2 — > move and further sim- 
plify the above expression as: 



,s,s 



^f^(-l) S<n ^ C «^ i * <(/3) ^ ! * <(5 ' ) %fe' 



55' 



/ j ±J ='qsp,pa u 5l_ 
5' 



DG, 



ikl,p-y 
qsp.pcr 



(132) 



The inverse of the above 4 — > 1 move is also true, thus we 
proof the invariance of the partition function Eq. (129) 
for those triangulations with different number of vertices. 

We believe such a Majorana valued partition function 
can describe almost all the quantum double of spin Chen- 
Simons theory. 42 



F. Simplified fixed-point conditions 

We have seen that under certain extra conditions, such 
as the tetrahedral conditions (100) and the orthogonal- 
ity condition (102), the F-tensors and O-tensors do not 
depend on the ordering of vertices: F t J l ™^ a = F^^f , 
Of a ' Q/3 = Of ' Q/3 . In this case eqn. (92), eqn. (93), and 
eqn. (94) can be simplified. We find that 

j?ijm' ,a' P' ( j-iijm,af}\* c 

• / .^kln.xS \ r kln, x 5 I ~ °maP,m'a'0' i 
nx5 

• F kln'x& = Wh(3n N Jim' < 1 Or N k ml* < 1 Or 

N kjn . < 1 or N nU , < l,or 

Skjn" (X) + S n il*(S) = odd, 

N kjt * N Hn . N lts , 

pijm,a/3 pit7i,ipxpjkt,r]K 
/ j / j / j / j knt,nip lps,n-y lsq,5<fi 
t n—1 ip—1 k—1 

-1 N qmp , 

= (_)*,,„,. J2 Fg3f*F«Z%. (133) 
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• E E E ot^io^y = i, 

k,j a=l 0=1 

• Of'-"' 3 = 0, if N ik .j. < 1, or JV«. ifc < 1, 



or Sifc.j.(o;) + Si*j k (P) = odd, 

e^O*"' ^^ = J2 (_)'«»•&)[««»• (x)+-i«,» (•»)]>< 



prnjp,7]/3 pij'm.aX/^jn* ,~/x 
r kln,x>. r nlk,jS ^k* 



for all satisfying -ZV^j;* > 0. 

• ^A l (A i )* = l, . i ! = A** ^ 0, 

i 

L J 

• det(M /cilia) g) = dot 



(134) 




^ 0. (135) 



VII. CATEGORICAL FRAMEWORK 

To provide a conceptual understanding of our gener- 
alization of string-net model, we discuss briefly the cat- 
egorical picture which underlies earlier algebraic manip- 
ulations. Such a mathematical framework will provide 
more examples for our fermionic string-net Hamiltonians 
in Appendix B. 

A string-net or Levin- Wen Hamiltonian can be easily 
constructed using 6j-symbols from a unitary fusion cat- 
egory C. The elementary excitations of the model form a 
unitary modular tensor category (UMTC) £ , which turns 
out to be the quantum double Z(C) of the input category 
C. A priori, the output modular category £ is not nec- 
essarily related to the input category C. Therefore, it is 
conceivable that similar Hamiltonians can be constructed 
from some other algebraic data where the elementary ex- 
citations still form a UMTC, which is not necessarily a 
quantum double. This is explored in Ref. 6. In the pre- 
ceding sections, we generalize the string-net model by 
including fermionic degrees of freedom. 

The mathematical framework for such a generalization 
is the theory of enriched categories. 24 An enriched cate- 
gory is actually not a category, just like a quantum group 
is not a group. We will consider only special enriched cat- 
egories, which we call projective super fusion categories. 
The ordinary unitary fusion categories are enriched cat- 
egories over the category of Hilbert spaces, while projec- 
tive super fusion categories are enriched categories over 
the category of super Hilbert spaces up to projective even 
unitary transformations. 

To the physically inclined readers, the use of category 
theory in condensed matter physics seems to be unjusti- 
fiably abstract. We would argue that the abstractness of 
category theory is actually its virtue. Topological prop- 
erties of quantum systems are independent of the micro- 
scopic details and non-local. A framework to encode such 



properties is necessarily blind to microscopic specifics. 
Therefore, philosophically category theory could be ex- 
tremely relevant, as we believe. 



A. Projective super tensor category 

We use super vector spaces to accommodate fermionic 
states, and generalize the composition of linear trans- 
formations to one only up to overall phases — a possibil- 
ity allowed by quantum mechanics. The projective ten- 
sor category of vector spaces is the category of vector 
spaces and linear transformations composed up to over- 
all phases, and the category of super vector spaces is the 
tensor category of Z2-graded vector spaces and all even 
linear transformations. 

In the categorical language, a fusion category is a 
rigid finite linear category with a simple unit. Equiv- 
alently, it can be defined using 6j-symbols: an equiv- 
alence class of solutions of pentagons satisfying certain 
normalizations 26 . Fermionic 6j-symbols F % kin^\b m 
eqn. (30) with certain normalizations define a projective 
super fusion category if they satisfy fermionic pentagon 
equations eqn. (52). However, the setup used in this pa- 
per may only generate a subclass of projective super fu- 
sion category. 



B. Super tensor category from super quantum 
groups 

The trivial example of a super tensor category is the 
category of Z 2 -graded vector spaces and all linear trans- 
formations. More interesting examples of super tensor 
categories can be constructed from the representation 
theory of super quantum groups. 

Super quantum groups are deformations of Lie 
superalgebras. 43,44 Though a mathematical theory analo- 
gous to quantum group exists, the details have not been 
worked out enough for our application here. In litera- 
ture, the categorical formulation focuses on the invariant 
spaces of even entwincrs, while for our purpose, wc need 
to consider all entwincrs. In particular, we are not aware 
of work on Majorana valued Clebsch- Gordon coefficients, 
therefore, we will leave the details to future publications. 



VIII. SIMPLE SOLUTIONS OF THE 
FIXED-POINT CONDITIONS 



In this section, let us discuss some simple solutions of 
the fixed-point conditions (91, 92, 93) for the fixed-point 
gfLU transformations (N ijk , N? jk , Si , F^ff, 0^). 
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A. Unimportant phase factors in the solutions 



B. A Z3 bosonic fixed-point 



Formally, the solutions of the fixed-point conditions 
are not isolated. They are parameterized by several con- 
tinuous phase factors. In this section, we will discuss the 
origin of those phase factors. We will see that those dif- 
ferent phase factors do not correspond to different states 
of matter (it different equivalence classes of gfLU trans- 
formations) . So after removing those unimportant phase 
factors, the solutions of the fixed-point conditions are 
isolated (at least for the simple examples studied here). 

We notice that, apart from two normalization con- 
ditions, all of the fixed-point conditions are linear in 
O^. Thus if (F^f/,O% a0 ) is a solution, then 



(Kin,iT/> elVa °i,a ) is also a solution. However, the 



two phase factors, e"* 



± do not lead to different 



fixed-point wave functions, since they only affect the to- 
tal phase of the wave function and are unphysical. Thus 
the total phases of 0^ Q ^ and CH^ a/3 can be adjusted. 
We can use this degree of freedom to set, say, O j. > 0. 

Similarly the total phase of X^;™'"^ is also unphys- 
ical and can be adjusted. We have used this degree of 
freedom to obtain eqn. (42) and eqn. (54) which do not 
contain the over all phase factor. 

The above two phase factors are unphysical. How- 
ever, the fixed-point solutions may also contain phase 
factors that do correspond to different fixed-point wave 
functions. For example, the local unitary transformation 
e i9i M io does jjQt a ff ec t the fusion rule Nijh, where M; is 
the number of edges with |?o)-state and |ZQ)-state. Such a 
local unitary transformation changes {F^^f^^ 1 ,0^ a ^) 
and generates a continuous family of the fixed-point wave 
functions parameterized by 9i . Those wave functions 
are related by local unitary transformations that contin- 
uously connect to identity. Thus, those fixed-point wave 
functions all belong to the same phase. 



Similarly, we can consider the following local unitary 
transformation \a) — > pU^ ^ that acts on each 
vertex with states |io), |jo)> l^o) on the three edges con- 
necting to the vertex. Such a local unitary transforma- 
tion also does not affect the fusion rule iVy^. The new 

local unitary transformation changes {F^i^^C -i^i^a ) 
and generates a continuous family of the fixed-point wave 
functions parameterized by the unitary matrix jj^ ok °\ 
Again, those fixed-point wave functions all belong to the 
same phase. 

In the following, we will study some simple solutions of 
the fixed-point conditions. We find that, for those exam- 
ples, the solutions have no additional continuous param- 
eters apart from the phase factors discussed above. This 
suggests that the solutions of the fixed-point conditions 
correspond to isolated zero-temperature phases. 



Let us consider a system where there are three states 
|0), |1), and |2) on each edge of the graph. We choose 

0* =0, 1* = 2, 2* = 1. (136) 

The simplest bosonic fusion rule that satisfies (91) is 

^000 = "ill = ^222 = -^120 = ^201 = ^012 
= ^210 = -X)2i = Ni$2 = 1, 
other Nijk = 0, 



N. 



f 



ijk 



0. 



(137) 



Since N^k < 1, there is no state on the vertices. So 
the indices a,/3, ... labeling the states on a vertex can 
be suppressed. The above fusion rule corresponds to the 
fusion rule for the string-net state that describes the Z 3 
gauge theory. 19 So we will call the corresponding graphic 
state Z 3 bosonic state (or simply Z 3 state), since 

N ljk = 

0. 

After fixing the unimportant phase factors, we find 
that the fixed-point conditions have only one isolated so- 
lution: 
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X- 



xx 
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Xi2°X+) { 
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X 



F, 
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210,+ 
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X = i 
X = i 



= 1 



^ir.H X. = 1 



Q 22 



- o 02 - o 11 - o 20 

1 

~7s' 



(138) 



We note that F^™'^/* in eqn. (25) relates wave func- 
tions on two graphs. In the above, we have drawn the two 
related graphs next to the F tensor, where the first graph 
following the F corresponds to the graph on the left-hand 
side of eqn. (25) and the second graph corresponds to the 
graph on the right-hand side of eqn. (25). The doted line 
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corresponds to the |0) -state on the edge and the solid line 
corresponds to the |l)-state or the |2)-state on the edge, 
depending on the direction of the edge. 

Fkinb can be determined from Fj^'_ (see eqn. (44)). 
It turns out that F l k \™'^ = Fj£™'_ for the Z 3 state. 
Similarly, o{ k Z a ^ can be determined from 0|V*^ (see 



eqn. (67)) and we find O 



jk,af3 



o 



Also the phase 



factors are all equal to 1 for the Z3 state: 



1. 



(139) 



We have introduced A 1 to describe the fixed-point wave 
function. The Z3 fixed-point wave function is described 
by 



00 



Q- 



10 



Q- 



20 



01 



a- 



11 



21 



02 



12 



22 



A> 



A 1 = A 2 



1 

7T 
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The above solution {F^^^^ 1 ,0{ k ^ a ^) gives rise to a 
fixed-point state which is an equal weight superposition 
of all the graphic states that satisfy the fusion rule N^k = 
1 on every vertex, where (|i), |fc)) are the three states 
on the edges that connect to the vertex. If we view the 
state |1) and |2) on a edge as a segment of string on 
the edge pointing in opposite directions and state |0) as 
no string state, then the graphic states that satisfy the 
fusion rule N^k = 1 can be viewed as states of orientable 
loops of strings. The fixed-point state is an equal weight 
superposition of all the loops of strings. 

If we put fixed-point state on a torus (ie a planar 
graph with periodic boundary condition in both x- and 
y-directions) , there are will be 9 states. This can be seen 
by viewing the orientable strings as flux lines carrying 
a flux quantum. The fusion rule Nijk implies that the 
flux is conserved mod 3. So the total flux quanta flowing 
across a non-contractible circle of the torus can be 0,1,2. 
There are two non-contractible circles on a torus which 
leads to 9 states (see Fig. 12). 

The Z3 state has a topological order described by the 
Z 3 gauge theory at low energies, 19,45 or equivalently by 
a U(l) x U(l) Chern-Simons (CS) theory: 46 



1 



£ = ^K IJ ai ^aj^ x + 



(141) 



with K 



3 
3 



FIG. 12. (Color online) The 9 fixed-point states on a "simple 
torus" . (The inner square and the outer square are identified 
to form a torus.) Any torus formed by more complex graph 
can always be reduced to the "simple torus" by the F-moves 
and the O-moves. The dotted lines represent the |0}-state on 
the edges, and the oriented solid lines represent the |l)-state 
or |2)-state on the edges. The pair of integers mn at the lower 
right corner of each graph indicates the total flux quanta (mod 
3) flowing across the two non-contractible circles of the torus. 



The simplest fermionic fusion rule that satisfies (91) is 
Nqqq = Niu = N222 = N120 = N201 = N012 

= N210 = Nq21 = N102 = 1, 

other N ijk = 0. (143) 



Aff — Aft — 1 
JV 111 — JV 222 — 

other N[ jk - 



0. 
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The above N^k is the same fusion rule for the Z3 string- 
net state discussed above. However, N(- k ^ 0. So we will 
call the corresponding graphic state Z3 fermionic state. 

After fixing the unimportant phase factors, we find 
that the fixed-point conditions have only one isolated so- 
lution: 



F 



000,+ ' 



000,- 



= 1 



200,-..' 



= 1 



C. A Z3 fermionic fixed-point 

Now let us discuss a fermionic example. We again con- 
sider a system with three states |0), |1), and |2) on each 
edge and we also choose 



0* = 0, r 
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F % kinb can be determined from F^™^ (see eqn. (44)). 
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It turns out that F 1 ^^ 



r kln,- 



for the fermionic Z3 
Similarly, 0{ can be determined from Oj k _l a ^ 

jk,aj3 



Of' a/3 . Also the 



state. 

(see eqn. (67)) and we find 0\ K ^ 
phase factors are all equal to 1 for the fermionic Z3 state: 
e i0£ _ e i0 O i — e i0o2 — 1 Actually the above fermionic 
Z3 state also satisfies the tetrahedral conditions (100), 
the orthogonality condition (102) as well as the associa- 
tivity condition (119). Thus such a solution is also an 
example of our fermionic generalization for Turaev-Viro 
invariance on closed triangulated 3-manifolds. 

The above fermionic solution {F^^ff^ 1 , 0- t' Q/3 ) give 
rise to a fixed-point state which is a superposition of all 
the graphic states that satisfy fusion rule Nijk = 1 on 
every vertex. However, the amplitude of each allowed 
graphic state can be +1 or —1. Using a similar argument, 
we can show that there are 9 fixed-point states on a torus, 
indicating that there are 9 types of quasiparticles for the 
fermionic Z 3 state (see Fig. 12). 

The fermionic Z3 state may be described by a U(l) x 

U{\) CS theory (141) with K = ^ and m = 

0,1,2,3,4,5. When m =even, the U{\) x {7(1) CS the- 
ory describe a bosonic theory (such as the m = case 
discussed before). In the following, we like to argue that 
the fermionic Z3 state corresponds to m = 3. 

Let us consider a Dehn twist of the torus. Here we 
define the Dehn twist formaly by requiring it to map a 

/ 

fixed-point state ^^(i, j, fc, a,/3) = 




\ 



on a torus to another fixed-point state defined on a differ- 



ent graph j, k, a, ft) = ip£~ 



3 



with 



^fix~(^.7Aa,/3) = ^~{i,j,k,a,P). Using a F-iaovc 
we can deform the graph 



3 



to 




which 



leads to an unitary transformation T between the fixed- 
point states on torus. The unitary matrix T can be cal- 
culated as the following: 



a_8 
fix~ 



n j ' k\ J 



a/3 
V'fix" 



E ijk,a8,af3 xl 



E ijk,aB,aJ3_ xl 
fi*jl*, X 6, x^fix 

i,xs 



3 
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Dip 



i' j' k' ,a' P' 



ijk,aj3, a/3 f oc 0' 

fix _ 



V 




If we choose the vertex label as a = a' = 1 and /3 
= 2, we obtain 



(\ 



a/3 
fix _ 



(\ 
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fix - 
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fix 
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This allows us to obtain: 
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The K = U ^ C/(l) x J7(l) CS theory has 9 types of 
quasiparticles with statistical angles 

(e ie *) = (1, 1, 1, 1, -1, e 1 ^ 3 , e i4 */ 3 , e i2T / 3 , e i5 ^ 3 ). 

We note that K — ^ gj CS theory is equivalent to the 

K = ^ CS theory. However, for the K = f jj ^ 

CS theory 

(e ie = (1, 1, 1, -1,-1, e iw / 3 , e i4 "/ 3 , e i2 " /3 , e i5w/3 ). 

In fact quasiparticle statistical angles e lBi always have 
± ambiguity in fermionic CS theory, since the fermionic 
electron is reguared as "identity". However, e 2l0i is well 
defined. We see that 

(e 2ie = (1, 1, 1, 1, 1, e i2 ^ 3 , e i2w / 3 , e i4w / 3 , e i4T / 3 ) 



fix 



^112,+ "rfix 



match the eigenvalues of T 2 exactly. 

For the K = ^ ^ U(l) x [7(1) CS theory, 

(e 2i9 = (1, 1, 1, 1, 1, e i27r / 3 , e i2 */ 3 , e i47r / 3 , e i47r / 3 ) 



lix 



= ( F 120.+ )V'fix 



where F's in ( ) are equal to —1 and other F's are equal 
to +1. We find that the 9 by 9 T-matrix for our fermionic 
Z 3 state is 



which also the eigenvalues of T 2 . However, the fermionic 
Z3 state cannot be the bosonic Z3 state described by the 
3^ 



K = I g [7(1) x U(l) CS theory. This is because A 1 

for the bosonic Z 3 state are all positive while A 1 and A 2 
for the fermionic Z3 state are negative. 
'0 3^ 



For the K — 



fj 17(1) x 17(1) CS theory, 
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For the bosonic states, the eigenvalues of T corre- 
spond to the statistical angles e l9i of the correspond- 
ing quasiparticle. 8 ' 23 For the fermionic states, we have a 
weaker result: 8 the eigenvalues of T 2 correspond to the 
square of the statistical angles, e 2l0i . 

The eigenvalues of the matrix T are 



(1, 1, 1, 1, 1, e l27T '\ e 12 ^ 3 , e 14 ^ 3 , e 14 ^ 3 ) 



,1271-/3 p i2rr/3 p i4Tr/3 A i47r/3 

and the eigenvalues of the matrix T 2 are 



(l,l,l,l,l,e i2 ^ 3 ,e i2 ^ 3 ,e i4w / 3 ,e i4 ^ 3 ). 



IOtt ■ 16ti 



(e zitt *) =(l,l,l,e 1 T,e 1 T 1 e 1 T ) e'T ) e'T ) e' 
For the K = ( ° jj ) U(l) x 17(1) CS theory 



14-7T ■ 14tT 



(e^ lW *) = (1, 1, 1, e'T, e'T, e 'T, e'r, e'T, e 'T) 



For the l\ = [ ° 3 A ) U(l) x U(l) CS theory, 



(e 2ie *) = (l,l,l,e i t,e i f ) e 1 T ) e i T ) e 1± r ) e i T) 
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(e alfl ') = (l,l,l,e 
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17(1) x U(l) CS theory, 



,2* J2JT jSjr jSjr ; 14k ■ liw 

» , e f , e 9 , e 9 , e 9 , e 9 ). 



They are all different from the eigenspectrum of T 2 . 
Therefore, our fermionic Z3 state is described by the 



K = 



3 
3 3 



3 
-3 



[7(1) x [7(1) CS theory. 
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D. A class of "trivial" fermionic fixed points 

There is a class of fermionic solutions that can be con- 
structed from the bosonic solutions. For every bosonic so- 
lution (N ijk ,N{ jk = Q,F^ff,Ot^,A*) that satisfy 
the conditions (91, 92, 93, 94), we can obtain a fermionic 
solution (N- •! N f — N--u F ijra ' a ^° nJ k ' a l 3 Ai\ 

SOlUllOIl (I^ijk, l^ijk — ^tjkj-^kln^X^ ' u i,a i A ) 

that satisfy the same set of conditions. Here 
(^^O^i*)are given by 

r kln,f\,b ~ auJr kln,j\,b 

(>!"■"■ <■<>::;■' # = a\ (iso) 

For example, when N = and Nq 00 = -Wooo = 1j we 
obtain the following solution 

£,000,ll,a 7 AOO.ll aO i fiKi\ 

F 000,ll,6 = flfe ' °0,a = a > A =1. (151) 

Such a solution give rise to a fermionic fixed-point state 
which has one fermion on each vertex and no other dy- 
namical degrees of freedom on the edges and vertices. 
This is a trivial fermionic state. We can view such a 
state as a fermionic "direct-product" state. 

For every bosonic fixed-point state described by 
(N ijk ,N{ jk = Q,F^ff,Oi^,A% we can combined 
it with the above fermionic "direct-product" state. This 
leads to the corresponding fermionic fixed-point state 
described by {N ijk ,N( jk = N ijk ,F^ff ,6(^,1% 
So, those fermionic fixed-point states are trivial in some 
sense. 



IX. SUMMARY 

Using string-net condensations and LU transforma- 
tions, we have obtained a systematic understanding of 
a large class of topological orders in bosonic systems. 6,19 
An interacting fermionic system is a non-local bosonic 
system. So classifying topological orders in fermion sys- 
tems appears to be a very difficult problem. 

In this paper, we introduced fLU and gfLU transfor- 
mations, which allow us to develop a general theory for a 
large class of fermionic topological orders. We find that a 
large class of fermionic topological orders can be classified 
by the data (N ijk ,Nf. k ,l$g;$f,0$ af> ) that satisfy a 
set of non-linear algebraic equations (91), (92), (93), and 
(94). Such a result generalizes the string-net result 6 ' 19 to 
fermionic cases. If we impose more conditions such as the 
tetrahedral conditions (100), the orthogonality condition 
(102) as well as the associativity condition (119) into our 
theory, we can further generalize the famous Turaev-Viro 
invariance to its fermionic version. We believe such a 
generalization will enable us to describe all the quantum 
double of fermionic CS theory. We hope our approach is 
a starting point for a general theory of topological orders 
in interacting fermion systems. 
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Appendix A: The fermionic structure of support 
space 

To understand the fermionic structure of the support 
space Va, let us first study the structure of pa- Let \</>i) 
be a basis of the Hilbert space of the region A and \4>i) 
be a basis of the Hilbert space of the region outside of A. 
\ip) can be expanded by <g> |<^): 

\^)=Y, c d<t>i)®\h)- (Ai) 

Then the matrix elements of pa is given by 

(pAh = Y^CaYCfl. (A2) 

i 

For a fermion system, the Hilbert space on a site, Vi has 
a structure: Vi = V® © , where states in V® have even 
numbers of fermions and states in V^ have odd numbers 
of fermions. The Hilbert space on the region A, Va, 
has a similar structure Va = Va ® V\, where states in 
VjI have even number of fermions and states in V\ have 
odd numbers of fermions. Let \<t>i, a ) be a basis of Vj£. 
Similarly, the Hilbert space on the region outside of A, 
Va, also has a structure Va — V^ (B V^. Let |^ a ) be a 
basis of Ujj . In this case, \ip) can be expanded as 

W = J2 c i,<xa,0\<f>i,«) ® \hp)- ( A3 ) 

the matrix elements of pa can now be expressed as 

Since the fermion number mod 2 is conserved, we may 
assume that \tp) contains even numbers of fermions. This 
means C i a .^ 7 — when a + 7 = 1 mod 2. Hence, we 
find that 

{pA)i,a;j,0 = 0, when a + (3 — 1 mod 2. (A5) 

Such a density matrix tells us that the support space Va 
has a structure Va = V^ ® Vj±, where V% has even num- 
bers of fermions and V\ has odd numbers of fermions. 
This means that U g contains only even numbers of 
fermionic operators (ie U g is a pseudo-local bosonic op- 
erator) . 

Appendix B: Ideal Hamiltonian for the fixed-point 
states 

In the section V, we have constructed the 
fixed-point wave functions from the solutions 
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each vertex: 



FIG. 13. A honeycomb lattice. The vertices are labeled by 
v, hexagons by p, and links by I. 



(N ljk ,N^F^^f l0 Zr P ^) of the self con S1 s- 
tent conditions. In this section, we will show that 
those fixed-point wave functions on a honeycomb lattice 
(see Fig. 13) are exact gapped ground state of a local 
Hamiltonian 



H = 5> -Q v ) + £(1 -B p ) 



(Bl) 



Vg = V G ® {®vV qubit ) 



(B3) 



where (£> v goes over all vertices and V qu bit is the two di- 
mensional Hilbert space of an auxiliary qubit \I), I = 
0,1. So in the expanded Hilbert space Vq x , the states 
on a vertex v are labeled by |a) ® I — 0,1. Vg 
is embedded into Vq x in the following way: each vertex 
state | a) in Vg correspond to the following vertex state 
I a) <8> \sijk{a)) in Vq x , where we have assume that the 
states on the three links connecting to the vertex are \i), 
\j) , and |fc). So the new auxiliary qubit \I) on a vertex 
is completely determined by k, a) and does not rep- 
resent an independent degree of freedom. It just tracks 
if the vertex state is bosonic or fermionic. The |0) -state 
correspond to bosonic vertex states and the |l)-state cor- 
respond to fermionic vertex states. 

In the expanded Hilbert space, Q v in H acts on the 
states on the 3 links that connect to the vertex v and on 
the states |a) ® |7) on the vertex v: 



where sums over all vertices and sums over all 
hexagons. We will first discuss the unprojective case 
where 



--±,b=± 



1. 



(B2) 



In this case, all the equal-up-to-total-phase relation ~ 
become the exactly-equal relation =. 

The Hamiltonian H should act on the Hilbert space 
Vq formed by all the graph states. It turns out that it 
is more convenient to write down the Hamiltonian if we 
expand the Hilbert space by adding an auxiliary qubit to 



: \cyk 



® | J) if N ijk > 0, 1 = s l3k (a), 



Qv o \I) = otherwise. (B4) 

Clearly, Q v is a projector Q\ = Q v . The B p operator 
in H acts on the states on the 6 links and the 6 vertices 
of the hexagon p and on the 6 links that connect to the 
hexagon. However, B p operator will not alter the states 
on the 6 links that connect to the hexagon. 

The Bp operator is constructed from a combination of 
the F-moves and O-moves as shown below: 




27 



/ J \ a 



EST^ (n sa' , x P" \i -j: a ' sa ^" ^'^lA T b'sb,j"-y,j^y -t~c'sc,X"X, A" A j 

/ / / / / / IVg, xff" / y "6*j*fc'*,7" / 3',r^ y "c*fe*c'%A"7', X^f J ' d'l'd 1 * , I m" X' , fix 
sa'x/3" b'7"/3' c'A"7' d' fi" X' 




EV^ ///isa'.x/i" \f j-a'sa,0" P, fc's6,7"7, 7^7 wj'sc.A" A, A" A 

Z^ Z^ Z^ Z^ ^ a,xP" ' J ~b*j*b'*,j"l3',j"/3' J ~c*k*c'*,X"y,X"y' J ~d'l' 



d",/j."\',fj."\' X 



sa'x/5" b'Y'fS' c'A"7' d' fi" X' e' V /i' 



e*m* e f * ,is" ft' , v" fi' fi x 



■,sa',x/9" \f 7 -a'sa,/3"/3, /3^/3 b'sb,r'j,Yj_ j-c sc.X" X, A" A 




y^ y^ y^ y^ y^ yi (Pa,^ y^b*j*b'*,j''/3\j^^c*k*c'*.x"j',x^^d*^d^^"x',f^ x 



sa'x/3" b'7"/3' c'A"7' d' fi" X' e'is"fi' f'a"u' 



d'sd,^" fj,, ^ fi -r-e'se,v"u, u" v j 

^e*m*e'*,i/"^', i/"/^' /*7l*/'*,a'V, a'V fi x 



EV- V"^ \"^/Y1 sa '' x < 9 " s |t -p a ' sa ^" P^I^A b'sb,~f"~/,~{^ -pc'sc, X" X, A^A 

/ , Z^ Z^ Z^ Z^ Z^ a,x£l - 1 ^b*j»b'*,Y'0', YJV J 'c'k'c",X"Y, X"y> J ~d'l'd>»,n"X>, t j,"X> X 
sa'x/3" b'~(" P' c'A"7 ; d' fi v X' e'i/ v fj,' j'a"v' x' a ' 




e*m*e'* [l' , is" fx' f*n* f* ,a" is' , q?j/_ a*i* a'* -x' a ' 5 X /q ' 



EV" V"^ V"^ /^.sa'.x^'^t x- a ' sa '' 3 ''' 3 '0 J r &'sfc,7"7,7^7 ^c'scA'^A, A" a 

/ , / , / / / , / , Z^^ a,X/9" ' ^b'j'b'* <v"fl' ^"fl'^.fc.^. \>W \<<~iJ-j,,,j,* „,>v „='x'X 




sa'x/3" b'y"/3' c' X" j' d' ^ X' e'l/'^i' fails' x' a ' 



j"/3', ~/"0' - / c*fc*c'*,A"7 / , X , ~t l J d*l*d'*,p,"X', fi"X' 
( \ a- /' \ 

u, H , ^ t se.jy .J-.J sj ,ut at. at at *r\8 a,\ XflY 

e*m* e'* .is" fi' , is" fi' f*n* /'* ,a'V , a" a* j* a 7 * .x'a 7 , x' a/ , x'X 



sd,fi" fi, fi' 7 fi -,-e' se.v" is. is" v - I ~f'sf,a"at. a" a a,x X 




= E E E E E E ^;^t v ,/v- (<, i, 

a'a' b'/3' c'7' d'A' e>' /V 



where the Majorana-number valued matrix # is given by 

K>aa,bfl,c-y,dKenJv ■ , , x _ V V V V V V ('O SQ '' x/3 '' ^ 7r a ' sa ' /3 "' 3 ' ^ 

a'Q'^'/3',c'7',d'A',eV,/'iy' V 4 ' m ' ^ 2^1 2^1 2^1 2^1 Z^i Z^i^ a^x^ > b* j*b'* , 7 "/3', 7 

«X^" 7" A" /i" 1/" a" x' 

-c'sc.A" A, A^A d'sd^'V,^> -j-e'se^'V, i^V 7 -/ , g/,a"a, a" a 



(B5) 



b'sb./y"i, 7" 7 e sc.A'A, A" A a sa,^ M M ^-e' se,v" j/, i^V j-f s},a"a, a" a /n,s~a,xx 

c*k*c f * ,X"y' , A" 7' d*l*d f * A', A' e*m*e / *,i/"/^') i/"^' /*n*/'*,a"i/', a'V 7 a*i*a f * : x f ct\ x' a ' a ' , x'x 



(B6) 
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Now, let us choose a particular order of the vertices. We assume the two new created vertices (3" and \ satisfy 
X < /3" < a, P, 7, A, )JL, 2A When we do the six F-moves, we can choose = j3, Y = 7, A = A, y! = fi, r/ = v, 
g[ = a, and f3" — 7" = A" = //' = 1/' = a" = \' ■ For such a choice of ordering, the above becomes 

sxP" 7" y m" ^" a" x' 

*^c* fe*c'*,A 1 '7 / J 4- d*l*d'* A',4- 6*771*6'* ,z/"ju/,4- *' y*^,*^/* jQt"^^ a* j t*a / * ^'ck'^ ct'^— 

-EEEEEEE p*) 

sx/3" 7" A" m" w" a" x' 

/ n a'*s*,/3"x-i*/) s ajt* ifi) n s a'jb" {&') pa'sa, 0"/3,+ /)S bfcc » (7)/)S(,'fce'* (7') pfc's&,7"7,+ z)»cM* (A) fl s c > !d /» (A') pc'se.A" A,+ 
IVa*,- j r b*j*b",-f"P',+° y 7' r c»fc*c'*,A"7',+ P A 17 A' r d*l*d'* ,/t" A',+ X 

ns Jm[ . {p.)n s d' mc " (aO rpd' sd,n" h,+ n 8 ™/' M ,)«<,'„/'* ( v ) E,e'se,i/"i/,+ /j s /io* ( Q )/j s /'.a'* ("') t-./'s/,q"q,+ n s*a : x'x 
% %' - r e*m*e'*,i/>', + t7 H <V r f* n* f* ,a" is' , + U °L ° a' r a* i* a" , X 'a' , + U a' - 



If we replace the two Majorana numbers 6* e and by 
the same Majorana fermion operator 77 q, on vertex a, 

#0 and by 77/3, etc , 2? will become the operator _B p 
that we are looking for. Here, the Majorana fermion 
operator r) v acts on the auxiliary qubit \I) on the vertex 
v: fj v \I) = |1 — I}. Or more precisely, after giving all the 
vertices an order, fj v is given by 



and 



n 



Vv — 11 u v' 



(B8) 



where a v) I = x,y, z are the Pauli matrices acting on the 
auxiliary qubit on the u-vertex. 

Now, we can express B p as a sum of products of op- 
erators which act on the 12 links and the 6 vertices of a 



(framed) hexagon 



B P = £ (B9) 

s,x/3"7" A"p"i/"q"x' 

f)S,P"Xps,P" p s .7" p s iA" jps,fj," ps,u" ps,a" As,x'x 
w l r /3;7" r 7;A" r X^" E;a" a;x 2 ' 

where the Fp?L„ operator, acting on the vertex (3 and its 




three legs, is given by 



F 



Epa'sa,/3"P,+ 
r b*j*b'*,-y"/3', 

a' abb' j (5 (3' 




(BIO) 



V/3 



The operators 0\'^ x and O s ^ x x act on the a- j3 link. 
They are given by 



01 



a'*s*,/3"x\ 



(Bll) 



2 s ' x x 



= E<- x,> 



(B12) 



The operators £? p and Fp-y» act on the expanded 
Hilbert space V@ x with auxiliary qubits. We can reduce 
them to act on the original Hilbert space Vg if we redefine 
the Majorana fermion operator r) v to be 



where 



E(- 

ijka 



T)('T 



(B13) 



(B14) 



In the following, we like to show that the matrix 
6 is a projection matrix satisfying B 2 = B, which 
means that B p is a projector: B 2 = B p . First 

tjaa,b/3 X7 .dX.ea, fu / ■ - ; ; \ 1 j 

Wwy/W«'W , *' l,m ' n ' can be viewed as 
an operator that act within the support space Vijkimn 

of the fixed-point state in the region surrounded by 
the fixed boundary states i, j, k, I, m, and n. Remem- 
ber that the support space is defined s the following: 
Let *fi x (a, a, b, (3, c, 7, d, A, e, ^, /, v, i, j, k, I, m, n; T) 
be the fixed-point wave function on the graph 



where T represent the indices 



on the other part of the graph which is not 
shown. If we fix i,j,k,l,m,n,T the wave func- 
tion %,j,k,i,m,n;r(o>, <x,b,/3, c, 7, d, A, e, fi, f, v) = 
* fix (a, a, b, (3, c, 7, d, A, e, ^, /, z/, i, j, fc, m, n; T) can 
be viewed as a state in the support space Vijkimn- 
As we vary V but keep i,j,k,l,m,n fixed, the wave 
functions *ij,fc,;, m ,„ ; r(a, «, &, P,c,j,d,\,e, fi, f,v) will 
span the full support space Vijkimn- The relation (B5) 
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FIG. 14. Up is generated by four F-moves and one O-move, 
which turn a hexagon in to a tree graph. (Up)* is generated 
by one O-move and four F-moves, which run in the reversed 
direction. 



implies that the operator represented by the matrix 

^o'a^^r^V'/tS/'^C*'^*'^^^ is the identit y 
operator in the support space Vijkimn- However, the 

above discussion did not tell us how B acts on a state 

that is not in the support space. 

To understand how B acts on a state that is not in 
the support space, let us consider the dimension Dijkimn 
of the support space Vijkimn which can be calculated by 

A „ /' A /• 



deforming the graph 




into 




If \m 1 1 

through a gfLU transformation U. Under the saturation 
assumption, Dijkimn is equal to the distinct labels in the 

A Ji 
tj 



graph 



(with fc, I, m, n fixed): 



D 



ijklmn ^ ] ^ijt* ^tku* ^ulv* N v 
tuv 



We note that the operator B is realized by a combi- 
nation of F-moves and O-moves (see eqn. (B5)) which 
realizes the following transformation 




However, we can perform the same transformation via 
other combination of F-moves and O-moves. The self 



consistent conditions satisfied by the i^-tensor and ex- 
tensor ensure that all those different ways to transform 
between the two states lead to the same B matrix. In par- 
ticular, we can have the following two paths that connect 

m o ; 

m — p / 




FIG. 15. A quantum state on a vertex-labeled graph G. The 
four vertices are labeled by a, /3, \, and 8 and the six edges 
are labeled by i, j, k, I, m, and n. For example, a = 1, 
P = 2, x = 3, 6 = 4, i = 5, i = 6, k = 7, 1 = 8, m = 9, 
n — 10 correspond to one vertex-edge labeling scheme. States 
on each edge are labeled by I = 0, ...,N. If the three edges 
of a vertex are in the states i, j, and k respectively, then 
the vertex has Nijk states. The states on the Q-vertex are 
labeled by a = 1, Nijk- Note the orientation of the edges 
are point towards to vertex. Also note that i — > j — > k runs 
anti-clockwise. 



the two states (see Fig. 14): 





(B15) where U is a gfLU transformation. We find that 



B = W p CUp 



(B16) 



where C, acting on ^ 




is a dimension 



Dijkimn x D 'ijklmn identity matrix. Also Up , containing 
only one O-move (see Fig. 14), has a form Up = U1PU2 
where U\p. are unitary matrices and V is a projection 
matrix. So the rank of B is equal or less than Dijki mn - 
Since it is the identity in the -Dijfc; m „-dimensional space 
Vijkimn, the matrix B is a hermitian projection matrix 
onto the space Vij k i mn - 
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a" a" b"/3" c"7" d"A" e"/^" /"i/" 

= ^&$%$Mj.AhSi *. l, m, n) (B17) 



In the above calculation of the Z?, we first insert a bub- 
ble on the a-link. We may also calculate B by first in- 
serting a bubble on other lines. All those different cal- 
culations will lead to the same B matrix, as discussed 
above. 

We note that B Pl B P2 and B P2 B Pi are generated by 
different combinations of F-moves and O-moves. Since 
the two combinations transform between the same pair 
of states, they give rise to the same relation between the 
two states. Therefore B Pi and B P2 commute 

B Pl B P 2 = B P2 B Vl- ( B18 ) 

We see that the corresponding Hamiltonian if is a sum 
of commuting projectors and is exactly soluble. 

Now let us consider the projective cases where 
(e 10F , e 1001 , e l6 ° 2 ) are not all equal to one. Following 
the above discussion, we can show that B p ~ (Up^Up. 
By adding a proper phase in B pi we can make B p a 
hermitian projector B p = B p = B p . However, the B p 
operators on neighboring hexagons may not commute 

B Pi B P 2 = e i8 ^B P2 B Pi . (B19) 

So if the phase factors (e ie ^, e iSo1 , e ie ° 2 ) are not all 
equal to one (ie for the projective case), the correspond- 
ing fixed-point state is a ground state of a sum of projec- 
tors, which may not commute. Note that the projectors 
are unfrustrated and the ground state energy is zero. 

Appendix C: Putting fermions on both vertices and 
edges 

In the main body of the paper, we put fermions only 
on the vertices. Actually, the whole formulation for the 
fixed point wave functions can be generalized to systems 
with fermions on both vertices and edges. In this section, 
we will describe such an generalization. 

1. Quantum state on a graph 

Here we also label the edges by i, j, etc (see Fig. 15). 
Now among the N + 1 states on an edge, N* of them 
are fermionic. We introduce s,; to indicate which state is 
fermionic: 

Si = 1 — > |i) is fermionic, 

Si = — > \i) is bosonic. (CI) 



Despite the similar look between i index and i index, 
and between a index and a index, the two sets of indices 
are very different. For example, a index labels the ver- 
tices while a index labels the states on a vertex. In this 
paper, we very often use a to label states on vertex a, 
and i to label states on edge i. 

In the string-net approach, we have assumed that the 
above graphic states on two graphs are the same if the 
two graphs have same topology. However, since different 
vertices and edges are really distinct, a generic graph 
state does not have such an topological invariance. Here 
we will consider vertex-edge-labeled graphs (ve-graphs) 
where each vertex and each edge are assigned an index 
a and i. Two ve-graphic is said to topologically the same 
if one graph can be continuously deformed into the other 
in such a way that vertex-edge labeling of the two graphs 
matches. We will consider the graph states that depend 
only on the topology of the ve-graphs. 



2. The first type of wave function renormalization 

The first type of renormalization does not change the 
degrees of freedom and corresponds to a local unitary 
transformation. It corresponds to locally deform the ve- 

i j k i j „ k 

graph 3f to \&n . We note that the support space 




should have the same 



number of fermionic states. Thus eqn. (16) can be split- 
ted as 

N 

m=0 
N 

+ E (N^N[ ml , + N( m ,N b kml ,)s m 

m=0 
N 

= 2~2(K n * N i-m + N£ jn .Nf.J(l - s n ), 

N 

+ >Z(Kn*rtni + <-„*^««>«- (C2) 
n=0 
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JV 



E ( N jim*N b km i* + N f ]im ,Nl ml ,)s m 

-0 

JV 

E ( N Jim^l nl , + Nl m ,N b kml ,){l - s m ) 

m.=0 

= E«n*^™ + N kjn* N Li>n, 
n=0 
JV 

+ E«'"*^m + <n»^m)(l " *»)■ (C3) 
n=0 

We express the above unitary transformation in terms 



m=0 



JV 



of the tensor F^™^' 11 , where a, & = ±, i,j,k,... 
0,...,N, and a = 1, iVyfe, etc : 



JV JVfcj,i* JVni,* 

kjki,r(a,/3,m) ~ E E E - F fei™x*,xF^' fei ' r ^' 5 ' n ) 

n=0 x=l 5=1 

(C4) 



or graphically as 



fix 



(C5) 



where the vertices carrying the states labeled by 
(a, (3, Xi 8) are labeled by ( a, (3, x, §) and the edges car- 
rying the states labeled by k, I, m, n) are labeled by 
(ij 3 1 ki L W) n) (see Fig. 15). 

In general, we should write the F-tensor as 
Fkln,x$, x^s~k~l~n~ w hi c h depends the vertex and edge la- 
bels a, /5, x> <L ii ii jfe: L Z/lj Z/d- Since the F-move only 
change the edge m to the edge n. We can always set 
m = n. Since there is an unique canonical assignment 
of edge labels, we assume that the F-tensor does not de- 
pend on the edge labels i, j, k, I, m, n. As before, we 
also assume that the F-tensor only depends on the sign 
of /3 — a. It does not depend on how big is the difference 

I n I <n 7-i • T-,*i m ) a /3 \ a 

\p — a\. So the b -tensor can be written as b kln %s ^— . 

Since some values of a, 0, ... indices correspond 

to fermionic states, the sign of wave function depend on 
how those fermionic states are ordered. In (C5), the 

(w A (\ J ^\ 

wave functions ips* \JK and V'flx \§*n are 
obtained by assuming the fermionic states are ordered in 



a particular way: 

V'fix 



fix 



1 Ki) 



sn r j x i 

EV'ffx*"" 1 - 2 '" ( J IxSnVlVZ-) xSnr, 1% ... 



(C6) 



where sum over the indices on the vertices and on the 
edges and rji are indices on other vertices. 

Thus eqn. (C5) should be more properly written as 



a /3 m.. 



1 „J. 



JV N kjn * JV„ i( * 

V - ^ ijm,aP, a/3 x§.n. 

2 / 2 / t kln, X S,xl %x 
n=0 x=l <5=1 



\^ J , (C7) 



where the superscripts a (3 m... and x^n... describing 
the order of fermionic states are added in the wave func- 
tion. 

Since the sign of the wave function depend on the or- 
dering of fermionic states, the F-tensor may also depend 
on the ordering. In this paper, we choose a particular 
ordering of fermionic states to define the F-tensor as de- 
scribed by a/3... and x<5... in eqn. (C7). In such a canon- 
ical ordering, we create a fermion on the f3- vertex before 
we create a fermion on the a- vertex. Similarly, we create 
a fermion on the 5-vertex before we create a fermion on 
the ^-vertex. 

After introducing one Majorana numbers 9 a on each 
vertex a and 9 i on each edge i, we can rewrite (C7) in 
a form that will be valid for any ordering of fermionic 
states on vertices and edges. We introduce the following 
wave function with Majorana numbers: 





(C8) 

where the order of the Majorana numbers 9 a 9 p9 m ... is 
tied to the order a /3 m.. in the superscript that describes 
the order of the fermionic states. We see that, by con- 
struction, the sign of ^gx ( I docs not depend on 
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the order of the fermionic states, and this is why the Ma- 
jorana wave function ^'^''j does not carry the 
superscript a (3 m... 

We like to mention that (9^,9 p) an d (9 x ,9g) are 
treated as different Majorana number even when, for ex- 
ample, a and x take the same value. This is because 
a and x label different vertices regardless if a and x 
have the same value or not. So a more accurate notation 
should be 



fix 




(C9) 



where 9 a and 9 X are different even when a = %. But in 
this paper, we will drop the ~ and hope that it will not 
cause any confusions. 

Let us introduce the F-tensor with Majorana numbers: 



(CIO) 



• WflSfcmi* iP)as m as na s n u* {S)a s kjn* (x) „ijm,aP, a/3 



We can rewrite (C7) as 



3. The second type of wave function 
renormalization 



The second type renormalization can now be written 



afii^ip j k.. 




Q3k,aP i t... ( i 



(C12) 



where i a and ip label the edges that connect to a and 
/3 vertices respectively, and the ordering of the fermionic 
states on the vertices and the edges is described by 
af3i a ipjk.... Note that the ordering of j and k is fixed 
by making ip — > j — > k to run anti-clockwise. 

The wave functions in eqn. (C12) is defined with re- 
spect to the ordering of the fermionic states described by 
a/3i a ipjk.... Let us introduce 



' .uM =\e s r^ {a) 9 s f^ (fi) 9 s ^9%x 




(C13) 



and 



O 



jk,a/3 _ a Sik*i* ( a )n s i*jk(P) nsi nsi nSj a s ka s ir) jk,a 



(C14) 



We can rewrite (C12) as 

* fix ( - of k ^ Rx (cis) 

which is valid for any ordering of the fermionic states. 



N N kjn * jv„ 4I , / i j k\ 

-EE E^Xf*- ■ (en) 

n=0 x=l 4=1 \ 1 / 

Such an expression is valid for any ordering of the fermion 
states. 



4. Fixed-point conditions on the F-tensor and the 
O-tesnor 

We can derive the fixed-point conditions on the 
F-tensor and the O-tesnor in the similar as de- 
scribed before. The fixed-point conditions are a 
set of non-linear equations whose variables are Nijk, 

N I 3 k, s ^ F kln,l\b^ °i k a af3 > and the seven P hases 
(e ie F | a =±,6=±, e i0oi ,e ie ° 2 . Here we just list the fixed- 
point conditions below 
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N N 



• Si = Si*, • N k ji = Ni*j* k *, • N kji = N i»j*k*' • Nji m *N km i* = N kjn *Ni* ni , 

m—0 n—0 

N N 

. ( N jim* N kml* + Nj im *N[ ml ,)s m + (N b oim ,N f kml , + Nf im ,N b kml *)(l - s m ) 

rii—O m—0 
N N 



^ ijm', a 'P',aP ( ijm xa p,aP^_ . „i iro , a /3, a x* _ r _N S;Hm , (a)« w . (/3) F jkn, X 5,i 

Z^ 1 kln, X S,xS \ r kln, X S,xS ) ~ v rnal3 ,m' a> 0' , m \ r kln, X SM I ~\ I r l'i*m*,, 

r kln, X S,+ y P kln, X 5,+ ' ~ e 2-*f kln, X 8,- \* kln,xS,- > ' 



kln, X 5,b 



• Kinxfb" = when ^J™* < 1 OT ^kmi' < 1 OT N kjn* < 1 or Nna- < I, or 



Sjim* (a) + s fem ;» {fi) + s kjn , (x) + s m i* (8) + s m + s n = odd, 

E(_)st[sti„* (f) + sits* («)] Tpijm,aP,+ pitn,tpx,+ pjkt,riK,+ 
V / knt,Tiip,— lps,K~f,— lsq,84>,— 



(_y jtm *(a)s lkq ,(S)+s m [s lnp ,( X )+s n ]+s q [s slp ,(~f)+s 3 +s m ] p,mkn,Px,+ T?ijm,ote,+ (017) 



. Qjk,afi ^Qjk,aPy _ ^ m Qjk,aP _ ^_y ik . (a)s jki . (0) + s i +s j s k Qk*j\Pa 

k./j a=l /3=1 

• i,t' Q ' 3 = °> if N ik-f < 1, or iVj^-fc < 1, or Sife.j-»(a) + Si* jk (P) + s< + Sj + s fe = odd, 

- „i0 O 2r) fc * J *'' 3Q! — ^ (ct)s jfci « (/3)+SjS fe \ " /■_\«4 1 ,*i«(7)»J»4*(A)+a4«» pij*k,a0,- p j'jm,i//i,+ n jn,J* 

9 e "^1*,+ — k / v / r jim,vn,+ J7 i*i*n*,A7,- ly i,4- ' 

m/ii/;n7A 

• e ie ° ± CP- n l' ct ' n 5ip8p& = y~] ^_yk jn *{x)ls n] * k 'h)+s zm * j *(a)]^_y z [i+s j +s m ]+s n is k +s m ]+s J is k +s m +s n ] x 

(_).«,. («^+^ +s ™+ s "]F™^F^X A ' + °i"*-; 7X ' for a11 fc ' ^ satisfying iV fcii . > 0. (C18) 



Finding N ijk , N( jk , Si , F%™ffi a , and Oj^ that sat- (N ijk ,Nf jk ,s i ,F^^ a ,O J i ^ a0 ) give us a characteriza- 

isfy such a set of non-linear equations corresponds to tion of fermionic topological orders (and bosonic topolog- 

finding a fixed-point gfLU transformation that has a icaJ orders as a special case with N f^ = Q and s . = ). 
non-trivial fixed-point wave function. So the solutions J 
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